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INTRODUCTION 


The  purpose  of  this  study  is  to  evaluate  the  design  of  the  industrial 
Liquid  waste  treatment  plant  presently  under  prove  out  at  Hols  ton  Army 
Ammunition  Plant .  ,  The  initial  design  work  was  completed  in  1974  by  Clark, 
Dietz  &  Associates  Engineers.2’3  Based  upon  their  preliminary  design  and 
additional  pilot  plant  studies  conducted  by  Holston  Defense  Corporation  in 
1977,  a  final  concept  design  was  developed  in  1977. 

The  major  concern  (in  this  study)  was  to  validate  the  present  design  and 
to  establish  guidelines  for  development  of  design  criteria  for  future  plants 
requiring  the  same  types  of  processes.  This  was  to  be  accomplished  by 
evaluating  data  collected  during  previous  studies,  and  then  performing  com¬ 
puter  simulation  studies  of  these  processes.  The  computer  simulations  were  to 
provide  information  on  the  transient  response  of  the  processes  due  to  varia¬ 
tions  in  the  mass  and  hydraulic  loading. 

LITERATURE  REVIEW 


GENERAL 

The  industrial  liquid  waste  treatment  plant,  which  is  shown  in  Eigure  1, 
consists  of  anaerobic  denitrification  filters  (anoxic  filters),  aerobic  fixed 
film  reactors  (trickling  filters),  and  aerobic  suspended  growth  reactors 
(activated  sludge) .  Only  the  wastewater  from  Area  B  is  to  be  treated  in  the 
anoxic  filters,  with  Area  A  wastewater  to  enter  the  system  just  ahead  of 
aerobic  fixed  film  reactors  (trickling  filters).  The  wastewater  from  the 
activated  sludge  processes  is  clarified  and  filtered  through  dual  media 
filters,  with  this  effluent  receiving  reaeration  prior  to  discharge  Into  the 
Holston  River.  The  piping  for  this  plant  is  such  that  the  activated  sludge 
system  may  be  operated  ahead  of  the  trickling  filters  or  in  the  normal 
sequence  with  the  trickling  filters  preceding  the  activated  sludge  process. 
Sludge  wasted  from  the  activated  sludge  process  is  thickened  before  being 
aerobically  digested  and  mechanically  dewatered.  The  dried  sludge  is  disposed 
of  in  landfill. 

The  characteristics  of  the  wastewater  being  treated  are  shown  in 
Table  1.  One  of  the  major  concerns  was  the  abiLity  of  this  plant  to  ensure 
that  the  effluent  complied  with  NPDES  waste  load  allocation.  The  area  of 
particular  concern  was  the  total  nitrate-as-nitrogen  loading.  The  1980  NPDES 
permit  is  summarized  in  Appendix  D. 


TABLE  1.  SUMMARY  OF  DESIGN  WASTEWATER  CHARACTERISTICS 
HOLSTON  ARMY  AMMUNITION  PLANT 


Parameter 

Design  Value  (Full  Mobilization)3 

Average  Flow 

12.6  MGD 

Maximum  Flow 

16.3  MCD 

COD 

61,535  pounds/day 

bod5 

43,000  pounds/day 

nh3~n 

510  pounds/day 

Total  P 

70  pounds/day 

Total  N 

3,375  pounds/day 

Total  SS 

4,092  pounds/day 

a.  Constructed  for  50%  of  mobilization  rates 


(assume  straight  line). 


Review  of  Denitrification 


Denitrification  is  coupled  with  the  aerobic  oxidation  of  organic  material; 
the  nitrate  which  acts  as  electron  acceptor  must  be  the  limiting  substrate 
present  in  the  system.  In  other  words,  the  carbonaceous  organic  concentration 
or  the  concentration  of  methanol,  which  may  be  used  as  an  organic  carbon 
source,  must  exceed  the  nitrate  concentration  in  order  to  effectively  deni¬ 
trify  the  wastewater.  Although  the  literature  implies  that  denitrification 
systems  may  be  designed  such  that  the  two  substrates  will  be  rate  limiting,  it 
is  difficult  in  practice  to  accomplish  this.  If  the  Monod19  kinetics  were  to 
be  used,  there  would  be  a  specific  growth  function  for  each  of  the  substrates 
as  shown  below. 
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where  p,  Sn,  Sc  are  the  maximum  rate  constant,  nitrate  substrate  concentra¬ 
tion,  and  carbonaceous  substrate  concentration,  respectively.  The  and  Kc 
are  the  half  rate  constants,  i.e.,  the  half  rate  constant  is  the  substrate 
concentration  at  half  the  maximum  growth  rate,  expressed  in  units  of 
mass/volume . 


In  practice  however,  the  organic  substrate  would  be  in  such  concentration 
(Sc»Kr)  that  the  nitrate  would  in  fact  be  the  rate  limiting  substrate,  and 
the  rate-of-growth  equation  reverts  to  a  saturation  function  involving  only 
nitrate,  as  shown: 
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While  the  suspended  growth  system  for  denitrification  represents  a  more 
efficient  process,  Grady4  found  that,  due  to  the  filamentous  nature  of  the 
biomass,  the  suspended  growth  system  could  not  be  used.  Therefore,  it  was 


recommended  that  the  submerged  anaerobic  filters  be  used.  Grady4  proposed 
that  the  system  followed  a  first-order  kinetics  expression  and  the  design 
could  be  based  upon 
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where  N  is  the  effluent  nitrate  concentration  in  milligrams  per  liter,  N  is 
the  influent  nitrate  concentration  in  milligrams  per  liter,  Z  is  the  filter 
depth  in  feet  and  V  is  the  hydraulic  application  rate  in  cubic  feet  per  second 
per  square  foot.  The  coefficient  lambda  was  found  to  be  equal  to  1.45  while 
the  power  of  the  hydraulic  application  rate  n  was  found  to  be  equal  to  0.5. 

This  expression  implies  that  the  removal  of  nitrate  is  essentially  a  logarith¬ 
mic  function  of  the  height  of  the  column.  Hash,  Evans,  and  Simerly5  showed  in 
their  work,  however,  that  there  were  two  distinct  rates  of  reactions  occurring 
within  the  column.  Through  the  first  4  feet  of  the  column,  the  reaction  was 
fast;  it  then  proceeded  at  a  much  slower  rate  through  the  next  (upper)  10  feet 
of  the  column.  These  data  are  shown  in  Figure  2.  Their  report5  indicates 
that  performance  of  the  anoxic  towers  (submerged  anaerobic  filters)  was  purely 
a  function  of  application  rate  and  that  the  kinetics  for  the  system  could  not 
be  clearly  described.  Harremose20  found  that  fixed-film  reactors  did  not 
follow  classical  first-order  kinetics  but  rather  followed  a  series  of  first- 
order,  half-order,  and  zero-order  reactions  which  overall  appeared  similar  to 
Monod  kinetics.  He  found  that  there  were  two  distinct  break  points  within  the 
kinetic  relationships  in  which  the  substrate  transport  as  well  as  a  reaction 
resistance  became  limiting,  as  shown  in  Figure  3. 

Kornegay  and  Andrews21  found  that  the  reaction  rate  was  a  function  of  the 
degree  of  penetration  of  the  substrate  into  the  fiLm  and  that  the  shear  effect 
of  the  fluid  had  a  significant  effect  upon  the  rate  of  overall  reaction.  It 
is  clear  from  all  the  information  presently  available  that  there  was  a  dis¬ 
crepancy  between  the  data  collected  by  Grady4  versus  that  of  Hash  et  al.5 
These  discrepancies  in  data  may  well  have  been  due  to  the  aging  of  the  samples 
forwarded  to  Grady  and  the  formulation  of  the  wastewater  used  in  his  experi¬ 
mental  work.  Hash  et  al.5  used  water  directly  from  the  production  process, 
which  was  representative  of  the  wastewater  that  will  be  treated  at  Holston. 

REVIEW  OF  TRICKLING  FILTER  DESIGN 

The  pilot  plant  work  conducted  by  Grady4  and  Clark  and  Dietz3  shows  the 
activated  sludge  process  to  be  unstable  due  to  development  of  filamentous  * 

growth.  Grady4  recommended  that  a  trickling  filter  or  attached  growth  system 
be  used  for  treating  the  carbonaceous  portion  of  the  wastewater.  Again, 

Grady4  was  able  to  describe  the  removal  of  carbonaceous  BOD  with  first-order  „ 

kinetics,  relating  the  rate  of  removal  with  depth  of  the  column  and  hydraulic 
application  rate.  The  expression  developed  was 
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reaction  rate  per  unit  surface  area 

zero-order  reaction  constant  per  unit  surface  area 

bulk  liquid  concentration 

saturation  constant 

(KfL2/D) l'2 

biofilm  constant  (mass  transfer  coefficient) 
depth  of  biofilm 
diffusion  coefficient 


Figure  3.  Trickling  filter  kinetics,  according 
to  Harremose  (20)  . 
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where  Z  is  the  depth  of  the  filter  and  Sc  and  SCQ  are  the  effluent  and  influ¬ 
ent  total  biochemical  oxygen  demand  in  milligrams  per  liter.  Lambda  was  found 
to  be  0.1  and  n  was  found  to  be  0.5  for  application  rate  (V)  in  cubic  feet  per 
second  per  square  foot.  Based  upon  this  study  it  was  recommended  that  four 
biological  filters  57  feet  in  diameter  with  a  media  depth  of  26  feet,  be  con¬ 
structed.  It  was  further  recommended  that  the  hydraulic  application  rate  be 
0.8  gallons  per  minute  per  square  foot. 

Hash  et  al.5  again  found  some  deviation  from  results  presented  by  Grady; 
they  found  that  lambda  actually  had  an  average  value  of  0.05  as  compared  to 
the  0.1  recommended  by  Grady.  They  agreed  that  the  hydraulic  application  rate 
exponent  was  essentially  equivalent  to  one-half.  This  is  consistent  with 
Rich22  where  n  is  actually  an  effectiveness  factor  for  the  hydraulic 
application  rate. 

At  similar  application  rates,  Hash  et  al.1’  observed  a  30  percent  reduction 
in  COD  as  compared  with  the  90  percent  reduction  reported  by  Grady.4  It  is 
difficult  to  explain  why  such  a  wide  deviation  would  exist  except  that  the 
characteristics  of  the  feed  streams  must  have  been  different. 

REVIEW  OF  ACTIVATED  SLUDGE  DESIGN 

Hash  et  al.5  further  investigated  the  possibility  of  using  activated 
sludge  as  a  treatment  process.  Applied  to  raw,  high  strength  wastewater,  the 
activated  sludge  system  gave  good  organic  removal,  but  poor  sludge  settling 
characteristics.  However,  activated  sludge  operating  at  a  solids  retention 
time  (SRT)  of  8  days  provided  an  excellent  second  stage  treatment  of  the 
trickling  filter  effluent.  The  results  (Table  2)  show  that  the  process 
removed  76  percent  of  the  COD,  85  percent  of  the  BOD  and  90  percent  of  the  NH^ 
-N  using  a  6-hour  aeration  cycle.  The  sludge  volume  index  was  110  mL/gm, 
which  indicates  good  settling  characteristics  of  the  sludge.  An  aeration 
period  of  4  hours  gave  slightly  inferior  performance  overall,  and  reducing 
sludge  age  to  3  days  gave  serious  bulking  problems .  It  was  recommended  that 
the  activated  sludge  plant  be  designed  for  an  aeration  time  of  6  hours,  a  SRT 
of  8  days,  and  a  flow  of  12.5  million  gallons  per  day. 


TABLE  2.  PERFORMING  OF  ACTIVATED  SLUDGE  SYSTEM  (From  Hash  et  al.5) 
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DEVELOPMENT  OF  DYNAMIC  MODEL 


GENERAL 

The  approach  used  in  developing  these  models  was  to  use  equations  which 
could  be  related  to  actual  design  variables.  It  was  also  necessary  to 
represent  these  design  variables  in  the  form  of  differential  equations  or 
integral  equations  so  that  the  dynamics  of  the  processes  could  be  evaluated. 

DENITRIFICATION  ANOXIC  TOWERS 

A  block  diagram  for  the  development  of  the  unsteady  state  model  for  the 
anoxic  towers  is  shown  in  Figure  4.  This  block  diagram  shows  the  input  and 
output  functions  as  they  relate  to  the  differential  equations  contained  within 
the  block  diagram.  These  differential  equations  were  developed  from  the 
literature  and  reports  by  Hash  et  al.5  and  Grady.4  Analysis  of  the  data  of 
Hash  et  al.5  shows  that  increasing  the  hydraulic  loading  rate  did  not  lower 
the  percent  removal  of  nitrate.  These  results  imply  that  their5  system  may 
have  been  substrate-limited,  since  the  towers  did  not  show  a  breakthrough 
characteristic.  Based  upon  the  information  available,  a  rate  equation  was 
developed  using  the  lambda  that  Grady  reported,  but  dependence  on  the  depth  of 
the  column  was  considered  a  fractional  order  function.  This  is  consistent 
with  observations  that  Eckenfelder23  made  in  trickling  filter  operation.  The 
fractional  order  for  the  operational  loading  rate  was  obtained  as  shown  in 
Figure  2.  The  data  of  Hash  et  al.5  were  replotted  and  found  to  be  essentially 
the  same  as  reported  by  Grady4  and  by  Eckenfelder.23  While  the  rate  did 
change  and  was  significantly  different  in  the  first  4  feet  on  the  column,  the 
next  (upper)  6  feet  had  essentially  the  same  rate  of  removal.  For  this  reason 
an  attempt  was  made  to  fit  the  data  to  zero-order  reaction  rate  as  well  as 
half-order  reaction  rate.  This  analysis  has  not  been  completed  at  this 
time.  The  primary  concern  was  to  obtain  a  rate  contour  for  the  total  depth  of 
the  tower  by  analyzing  the  data  of  Hash  et  al.5  A  mean  value  of  1.20  for  the 
power  term  (R  in  equation  9  below)  of  depth  was  determined.  Using  this  term, 
the  Continuous  Simulation  Modeling  Program  ( CSMP)  equations  for  the  model  were 
written.  Modeiing  of  the  denitrification  system  was  done  on  a  macro  basis  and 
did  not  in  any  way  consider  characteristics  of  transport  within  the  biofiim. 
Transport  in  the  biofilm  is  being  looked  at  as  another  element  of  the  study, 
but  has  not  been  thoroughly  documented  at  this  point.  The  output  observed 
from  this  model  is  shown  in  Appendix  A.  The  general  structure  of  the  model  is 
as  follows: 

rate  of 
mass  N 
in 

[Fm<  -  Fm  «  ]  AXAY+[r»]  AXAYAZ  =  AXAYAZ  (6) 

N  Z  N  Z+AZ  *  9t 


rate  of  rate  of  rate  of 

mass  N  +  mass  N  =  mass  N  (5) 

out  reacting  accumulating 
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Figure  4.  Anoxic  tower  model. 


Dividing  Equation  6  by  (AX,  AY,  AZ),  the  differential  sections,  and  taking  the 
limit  as  AZ  +0,  the  following  equation  is  obtained 


or 


9Jj 

3Z 


+  r , 


3N 
3  t 


dFN  =  dN 

dZ  r  N  dt 


(7) 

(*) 


where  is  the  flux  term,  N  is  the  nitrate  concentration,  and  r^  is  the 
reaction  term. 


The  reaction  term  is  expressed  empirically  in  terms  of  Z 


dN 

CN  dt 


—  )  N 
(0/A)W 


(9) 


where  X  is  a  rate  coefficient,  R  is  an  effectiveness  factor,  and  w  is  a 
function  of  the  packing  used  in  the  column. 


The  function  Z  represents  the  column  depth,  while  Q  represents  flow  rate 
in  cubic  feet  per  second,  and  A  is  the  cross  sectional  area  of  the  column, 
which  gives  the  velocity  through  the  column. 


The  change  in  flux  may  be  expressed  as 


(10) 


afn  = 


mass 


length^  x  time 
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Therefore,  the  unsteady  state  equation  (Equation  8)  may  be  expressed  as 


VdN 

dt 


Q  (Nq-N)  +  rN  V 


(11) 


Then,  eliminating  the  volume  term  (V)  where  t  =  (V/Q) ,  the  equation  becomes 


dN 

dt 


-  i  ( »0-H> 


+  r 


N 


(12) 


This  differential  equation  describes  the  function  of  the  denitrification 
columns.  In  order  to  evaluate  the  dynamic  response  of  the  system,  the  program 
was  written  so  that  the  Influent  flow  and  nitrate  concentration  could  be 
varied  independently  with  time.  This  was  done  by  using  either  a  step  or 
impulse  function.  Figure  5  shows  how  these ■ functions  were  defined  within  the 
program. 


FUNCTION  TO  START  SQUARE  WAVE  PULSE 

TRIG  1  =IMPULS  (6. 0,6.0) 
y 


# 

t 

FUNCTION  FOR  MASS  LOADING 


SNO«SNIxPULSE  (6.0,TRIG1) 


Figure  5.  Transient  function  for  computer  simulations. 
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TRICKLING  FILTER  SYSTEM 


The  trickling  filter  model  was  developed  in  exactly  the  same  manner  as  the 
anoxic  towers  model.  Replacing  N,  the  nitrate  concentration  with  S,  the 
carbonaceous  substrate  concentration,  an  equation  for  the  trickling  filter  was 
expressed  as  follows: 


dS 

dt 

t  <  So‘S) 

+  rs 

(13) 

where  r 

'  *ZB  )  s 

(14) 

dt  1 

(Q/A)U 

The  parameters  Q,  A,  and  Z  represent  flow,  cross-sectional  area,  and  column 
height,  just  as  in  the  anoxic  tower.  The  terms  B  and  u  represent  an  effec¬ 
tiveness  factor  and  a  media  packing  factor.  As  in  the  case  of  anoxic  towers 
the  Q/ A  term  represents  a  velocity  and  the  u  represents  a  holdup  adjustment 
for  the  media.  Lambda  (A)  represents  a  rate  term. 

ACTIVATED  SLUDGE  SYSTEM 


The  activated  sludge  model  was  based  upon  Monod  kinetics.  Substrate  and 
organism  material  balances  were  written  over  the  entire  system,  as  follows: 


(SUBSTRATE  BALANCE) 

j  c  n  pn  X  S 

~  =  5  (  S  -S')  -  m  _ 

dt  Vo'  Y(S  +  K  ) 


m 


(ORGANISM  BALANCE) 


X  S 


#-! <VX>  +  s Vr-  V 

m 


(15) 


(16) 


where  SQ  and  S  are  influent  and  effluent  substrate  concentrations,  respec¬ 
tively,  XQ  and  X  are  influent  and  effluent  organism  concentrations,  respec¬ 
tively,  Y  is  the  mass  ratio  of  cells  formed  to  substrate  consumed,  and  p  ,  , 

and  Kd  are  constants. 

The  constants  used  in  the  model  were  those  recommended  for  design  in  ASCE 
Manual  of  Practice  Number  36. 24  Solid  retention  time  of  8  days  was  used  in 
order  to  tune  the  model  with  the  data  collected  at  Holston.  It  is  expected 
that  once  the  model  is  tuned,  the  SRT  will  be  varied  in  order  to  assess  the 
system's  ability  to  respond  to  transient  loadings. 

GENERAL  ASSESSMENT  OF  MODEL 

The  Continuous  Simulation  Modeling  Program  (CSMP)  is  capable  of  predicting 
the  effluent  discharge  from  the  Holston  AAP  industrial  waste  treatment 
plant.  The  model  is  based  upon  the  rational  design  equations  used  in  current 
practice  and  is  consistent  with  equations  used  to  design  the  wastewater  treat¬ 
ment  plant.  All  expressions  used  in  this  model  have  physical  significance;  no 
empirical  mathematical  expressions  were  used. 


* 
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The  model  is  also  written  in  BASIC  language  for  use  on  the  Apple^*  com¬ 
puter.  The  program  in  BASIC  appears  to  be  better  suited  to  the  process  than 
CSMP.  CSMP  did  not  adapt  well  in  modeling  the  trickling  filter,  due  to  the 
nature  of  the  equation  describing  the  trickling  filter.  The  CSMP  program  also 
has  the  disadvantage  of  requiring  a  special  compiler,  while  the  BASIC  program 
can  be  run  on  the  Applet  or  other  microcomputers. 

In  the  verification  phase  of  the  study,  the  mass  transfer  rates  for  the 
anoxic  tower  and  the  trickling  filter  will  be  studied. 


DISCUSSION  OF  RESULTS 


GENERAL  DESIGN  COMMENTS 

Review  of  the  proposed  design  indicates  that  the  modifications  recommended 
by  Hash  et  al.5  were  justified.  These  modifications  to  the  original  design 
certainly  provide  greater  flexibility  in  the  operation  of  the  wastewater 
treatment  plant  and  a  greater  degree  of  assurance  in  meeting  the  NPDES  cri¬ 
teria.  The  development  of  the  simulation  models,  which  are  discussed  below, 
shows  that  all  processes  in  the  system  are  compatible  and  that  the  plant 
should  function  well  under  wide  variation  in  loadings. 

In  the  verification  phase  the  mass  transfer  rates  for  both  the  anoxic 
tower  and  trickling  filter  need  to  be  further  refined.  An  attempt  will  be 
made  to  correlate  the  operational  data  with  the  height  and  number  of  theo¬ 
retical  units  or  a  function.  If  this  correlation  can  be  developed  it  will 
permit  the  sizing  of  a  column  for  a  specific  efficiency  and  capacity  very 
quickly  and  easily.  The  simulation  studies  which  are  discussed  below  were 
based  upon  the  most  extreme  conditions  that  might  be  encountered.  Results 
indicate  that  the  industrial  waste  plant  has  the  capacity  to  respond  both  to 
extreme  organic  and  hydraulic  surges.  While  the  simulations  do  not  show  that 
the  system  complies  with  all  the  NPDES  criteria,  there  is  a  trend  which 
implies  compliance.  When  data  are  available  to  tune  the  models,  more  accurate 
predictions  may  be  made. 

CONTINUOUS  SIMULATION  MODELING  PROGRAM 

The  initial  modeling  of  the  wastewater  treatment  plant  was  done  using  a 
Continuous  Simulation  Modeling  Program  (CSMP).  This  provided  an  easy  means  of 
developing  the  differential  equation  in  a  computer  format.  This  model  is 
shown  in  Appendix  A.  The  dynamic  section  of  the  program  contains  the  anoxic 
towers  followed  by  the  trickling  filters  and  activated  sludge  processes. 

The  CSMP  graphically  displays  the  inputs  and  output  of  each  process.  For 
example,  the  anoxic  towers  with  inputs  ranging  between  25  and  75  mg/L  of 
nitrate  nitrogen  showed  effluent  concentrations  of  3  mg/L  to  15  mg/L.  The 
transients  in  the  tower  effluent  lag  the  input  transient  by  approximately  30 


*  Use  of  trademarked  name  does  not  imply  endorsement  by  the  US  Army,  but  is 
used  only  to  assist  in  identification  of  a  specific  product. 
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minutes.  Sensitivity  analysis  of  the  parameters  needs  to  be  performed;  when 
actual  operational  data  are  available  this  will  be  accomplished. 

Simulation  of  the  trickling  filter  showed  that  effluent  from  the  trickling 
filters  varied  from  16  mg/L  to  35  mg/L.  The  influent  pulse  ranged  from  2.50  to 
350  mg/L  with  the  effluent  response  lagging  by  1  hour.  These  simulations 
indicated  an  efficiency  of  the  filters  ranging  between  90  and  94  percent.  » 

These  results  were  slightly  different  than  those  given  by  the  BASIC  Simulation 
Program,  but  are  in  agreement  with  the  findings  of  Hash  et  al .5 

The  simulated  effluent  concentrations  from  the  activated  sludge  process 
showed  a  concentration  of  carbonaceous  material  which  ranged  from  5  mg/L  to 
50  mg/L.  The  organism  concentration  was  approximately  100  mg/L,  which  implied 
that  the  system  was  dispersed  and  substrate  limited.  Until  actual  data  are 
available,  it  will  be  difficult  to  determine  if  this  is  a  representative  case 
of  actual  operating  conditions. 

BASIC  SIMULATION  MODELING  PROGRAM 

A  simulation  of  48  hours  of  operation  was  made  using  the  Basic  program, 
shown  in  Appendix  B.  The  initial  nitrate  concentration  of  25  mg/L  was 
periodically  pulsed  to  50  mg/L.  Flow  rates  were  also  pulsed  from  15.0  million 
gallons  per  day  (mgd)  to  19.5  mgd  to  simulate  hydraulic  disturbances.  The 
output  from  the  simulation  (Appendix  B)  showed  that  the  program  was  responsive 
to  the  changes  in  the  input  parameters,  with  effluent  nitrate  concentrations 
ranging  from  3  to  12  mg/L.  It  is  expected  that  with  proper  tuning  of  the 
model,  greater  sensitivity  may  be  obtained.  Verification  of  the  model  with 
actual  plant  data  is  essential  to  the  development  of  an  accurate  effectiveness 
ratio  and  rate  constant. 

Tine  same  simulation  study  was  used  to  evaluate  the  trickling  filter  por¬ 
tion  of  the  model.  A  carbonaceous  organic  concentration  of  25D  mg/L  was  used 
in  the  simulations.  Using  the  rate  constants  and  effectiveness  factor  previ¬ 
ously  described  for  the  trickling  filter,  the  simulations  predicted  between  87 
and  91  percent  reduction  in  the  carbonaceous  concentration.  Based  upon  the 
limited  data  available,  these  results  seem  appropriate.  Again,  verification 
of  the  systems  using  actual  plant  data  will  provide  a  better  assessment  of  the 
rate  constants  and  effectiveness  coefficients. 

The  activated  sludge  process,  which  follows  the  trickling  filter,  was 
simulated  using  substrate  concentrations  which  ranged  between  50  mg/L  and  150 
mg/L.  Hydraulic  transients  were  also  used  in  thLs  process  analysis.  Results 
of  the  simulation  showed  that  effluent  would  range  between  5  mg/L  and  77  mg/L  * 

using  the  rate  constants  selected. 

As  indicated  throughout  the  discussion  of  the  modeling  phase  of  this  t 

study,  the  future  use  of  actual  operational  data  to  refine  the  rate  constants 
will  make  the  models  more  predictive.  These  models  have,  however,  proven  that 
the  processes  selected  in  the  design  are  interactive  and  should  prove  to  be  an 
effective  treatment  system.  Certainly  these  processes  warrant  consideration 
in  treatment  of  any  wastewater  containing  these  contaminants. 
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CONCLUSIONS 


The  following  conclusions  may  be  made: 

1.  Design  of  the  Holston  Industrial  Waste  Treatment  Plant  is  consistent 

*  with  good  engineering  practice. 

2.  Simulation  studies  using  the  model  previously  described  indicates  that 
this  plant  would  meet  the  NPDES  permit  criteria. 

3.  Transient  analysis  of  the  treatment  plant  shows  that  the  processes  are 
capable  of  sustaining  any  expected  upsets  without  exceeding  permit  criteria. 
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****l>JPf-hKENi  iAL  ANALOG  Di  ti  1  I AL  SYSfEI'l**** 

***PROBLEM  INPUT  STATEMENTS*** 

**************************  HOLSTUN  AAP  ******************************- 
*****************  INDUSTRIAL  WASTE  TREATMENT  PLANT  ******************* 
**************CARNAHAN,  MARCH  1 933************************* 

INITIAL 

*  ANOXIC  TOWERS 

CONSTANT  A=6655.  0,  ZLAM=1.  45,  U=0.  5,  R=0.  33,  .  .  . 

Z=10.  0,  VOL=66550.  0 

CONSTANT  QIN2=7.  5,  GIN3=2.  25,  SNINF=25.  O,  .  .  . 

SN1=50.  0,  SNI=25.  O 

*  TRICKLING  FILTER 

CONSTANT  CA=  1 0000.  0,  CLAM=0.  71,  ZT=26.  O,  W=0.  5,  .  .  . 

B=.  33,  0V0L.=260000.  0 

CONSTANT  0102=1.  1,  QIC3=0.  25,  SC  I NF=250.  O,  .  .  . 

SC  1  =  1 00.  0,  SC  I  =250.  O 

*  ACTIVATED  SLUDGE 

CONSTANT  K=5.  O,  KS=60.  O,  KD=0.  OS,  YG=0.  60,  .  .  . 

AV0L=2S70000.  0,  QR=0.  25 

***************  ASSUMED  XC  SHOULD  BE  0.  0  ***************** 

I  NOON  XI  =  100.  0,  SC =50.  0,  XC=100.  0 

DYNAMIC 

Kl=K/24 
K2=KD/24 
V 1 = VOL 
V2=CV0L 
V 3= A VOL 

*  ANOXIC  TOWER 

TR I G 1  =  I MPULS < 2.  0,  6.  0) 

TR  IG2=IMPUl.  S<  2.  0,  4.  O) 

TR  I  G3=  IMPULS  <  2.  0,6.  0) 

TR  IG4=IMPULS(  2.  0,  3.  0) 

QIN1  =  (QIN3/ <24.  0*60))*10**6 
Q I NF =  <  0 1 N2 / <  24 .  0*60) >*10**6 
0 I C 1 = ( Q I €3/ ( 24. 0*60) >*10**6 
QICF  =  (QIC2/ <24.  0*60) >*10**6 

SN0=SN1*PULSE<2.  0,  TRIG1 )+SNINF 
00=0 1 N 1  *PULSE  <  3.  0,  TR  I G2 )  +0 1  NF 
SCO=SC 1 *PULSE < 2.  0, TR IG3 ) +SC I NF 
0T0=QIC1*PULSE<2.  0,  TRIG4)+0ICF+Q0 
0T0A=QT0+ <  QR*OTO ) 

X0= ( <  X I *GT 0 )  +  <  X*QR*QTO ) ) /QTOA 
SAF=  < <  QTO/CA ) **W ) 
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*******  *:*-s*r***SECOND  SECTION********************** 


SNF-- 1 NTGRL  (  SN 1 ,  DOTSN ) 

DOTSN=<  (G0*3.  021 ) /VI ) *< SNO-SNF )  +  < < -ZLAM* < Z**R ) >/  <  ( QO/A ) **U ) ) *Sh 
SCF=  i  NTGRL  (SCI,  DOTBC  ) 

DDT3C=(  (GT0*8.  021  ) /V2 ) * < SCO-SCF )  +  (  (  -CLAM*  (  ZT **B  )  )/SAF)*SCF 


DSD'I  =  <  <  (QT0A*3.  021  )  /V3 )  *(  SCF-S )  )  -  <  K1  *  <  3/  (  3  +  KS )  )  *X ) 

DXDT- ( <  <GTOA*B.  021 )/V3)*(X0~X> >  +  <  < K1 *YG > *< S/ ( S+KS ) )*X)-(K2*X> 
S= I NTGRL (SC, DSDT) 

X=I NTGRL (XC, DXDT ) 

NOSORT 

IF(X.  GT.  10000.  0)  X=10000.  0 
IF(SCF.  LT.  10.  0)  SCF=  1 0.  0 
IF  (  S.  LT.  5.  0)  S-5.  O 


TERMINAL 


METHOD  MILNE 

TIMER  PRDEL=0.  5,  0UTDEL=0.  5,  FINTIM=48.  O,  .  .  . 

DELT=0.  1 

PRINT  GO,  QTOA,  QTO 

LABEL  INFLUENT  NITRATE  CONCENTRATION 

PRTPLT  SNO 

LABEL  EFFLUENT  NITRATE  CONCENTRATION 

PRTPLT  SNF 

LABEL  T.  F.  INFLUENT  CONCENTRATION 

PRTPLT  SCO 

LABEL  T.  F.  EFFLUENT  CONCENTRATION 

PRTPLT  BCF 

LABEL.  A.  S.  EFFLUENT  CONCENTRATION 

PRTPLT  S 

LABEL  ORGANISM  CONCENTRATION 

PRTPLT  X 

END 

STOP 


OUTPUT  VARIABLE  SEQUENCE 


TRIG1 

SNO 

VI 

QINF 

TRIG2 

QIN1 

GO 

DOTSN 

SNF 

QICF 

TRIG4 

GIC1 

QTO 

SAF 

TRIG3 

SCO 

V2 

DOTSC 

BCF 

K 1 

V3 

S 

QTOA 

DSDT 

S 

K2 

XO 

DXDT 

X 

X 

SCF 

OUTPUTS 

INPUTS 

PAR AMS 

INTEGS 

:  +  MEM 

BLKS 

FORTRAN 

DATA 

CDS 

35(500) 

66 ( 1400) 

35(400) 

4  + 

0=  4  <  300 ) 

32(600) 

2 

8 
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CONSTANT  A=6655  0.  ZLAM=1.45>  U=0.  5,  R=0.  33, 

Z=10.  0,  V0L=66550.  0 

CONSTANT  GIN2=7.  5,  QIN3=2.  25,  SNINF=25.  O,  .  .  . 

SN1-50.  0,  SNI=25.  O 

CONSTANT  CA=  10000.  0,  CL.AM=0.  71,  ZT=26.  0,  W=0 

B“.  33,  C VOL.— 260000.  0 

CONSTANT  GIC2=1.  1,  QIC3=0.  25,  SC  INF=250.  0,  .  . 

SC  1  =  100.  0,  SC  I  =250.  O 

CONSTANT  K=5.  0,  KS=60.  O,  KD=0.  06,  YG=0.  60,  .  . 

AV0L=2B70000.  0,  GR=0.  25 
INCON  XI  =  100.  0,  SC=50.  0,  XC=100.  0 

METHOD  MILNE 

TIMER  PRDEL=0.  5,  0UTDEL=0.  5,  FINTIM=48.  O,  .  .  . 
DELT=0.  1 

PRINT  GO,  Q  TOA,  QTO 

LABEL  INFLUENT  NITRATE  CONCENTRATION 

PRTPLT  SNO 

LABEL  EFFLUENT  NITRATE  CONCENTRATION 

PRTPLT  SNF 

LABEL  T.  F.  INFLUENT  CONCENTRATION 

PRTPLT  SCO 

LABEL  T.  F.  EFFLUENT  CONCENTRATION 

PRTPLT  SCF 

LABEL  A.  S.  EFFLUENT  CONCENTRATION 

PRTPLT  S 

LABEL  ORGANISM  CONCENTRATION 

PRTPLT  X 

END 

TIMER  VARIABLES 
DEL.T  =  1 .  2500E— 0 1 
DELHI M=  4.  S000E-06 
F I NT I M=  4.  OOOOE  01 
PR  DEL  =  5.  OOOOE— 0 1 
OUTDFL=  5.  OOOOE —01 
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n  lit: 

00 

0  to  A 

010 

o.  ooooe- 

-01 

5  2003E 

03 

7  4603F 

03 

fi.  9722E 

03 

5.  ooooe  - 

-01 

0.  2 08 LIE 

03 

7.  4603E 

03 

0.  9782E 

03 

1.  OOOCE 

00 

5.  208 3E 

03 

7.  4653E 

03 

5.  9722E 

03 

1.  OOOOE 

00 

5.  2003E 

03 

7.  4603E 

03 

5.  9722E 

03 

2.  OOOOE 

00 

6.  770GE 

03 

9.  6354E 

03 

7.  7083E 

03 

2.  5000E 

00 

6.  770SE 

03 

9.  6354E 

03 

7.  7083E 

03 

3.  OOOOE 

00 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

3.  5000E 

00 

6.  7703E 

03 

9.  6354E 

03 

7.  7083E 

03 

4.  OOOOE 

00 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

4. OOOOE 

00 

6.  7703E 

03 

9.  4 1 84  E 

03 

7.  5347E 

03 

5.  OOOOE 

00 

6.  7703E 

03 

9.  6304E 

03 

7.  7083E 

03 

5.  OOOOE 

00 

5.  2083E 

03 

7.  6323E 

03 

6.  145SE 

03 

6.  OOOCE 

00 

6.  7708E 

03 

9.  6354E 

03 

7.  7083E 

03 

6.  OOOOE 

00 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

7.  OOOOE 

00 

6.  7708E 

03 

9.  6354E 

03 

7.  7083E 

03 

7.  OOOOE 

00 

6.  7703E 

03 

9.  4184E 

03 

7.  5347E 

03 

8.  OOOOE 

00 

6,  770SE 

03 

9.  6304E 

03 

7.  7083E 

03 

8.  OOOOE 

00 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

9.  OOOCE 

00 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

9.  OOOOE 

00 

5.  2083E 

03 

7.  6323E 

03 

6.  1408E 

03 

1 .  OOOOE 

01 

6.  7703E 

03 

9.  6304E 

03 

7.  70B3E 

03 

1.  OOOOE 

01 

6.  7708E 

03 

9.  4184E 

03 

7.  5347E 

03 

1.  1000E 

01 

6.  770SE 

03 

9.  6354E 

03 

7.  70S3E 

03 

1.  1000E 

01 

6  7708E 

03 

9.  6354E 

03 

7.  70B3E 

03 

1.  200CE 

01 

6.  7708E 

03 

9.  6304 E 

03 

7.  7083E 

03 

1.  2000E 

01 

6.  7708E 

03 

9.  6354E 

03 

7.  7083E 

03 

1.  3000E 

01 

6  770SE 

03 

9.  6354E 

03 

7.  70S3E 

03 

I . 3000E 

01 

0.  2083E 

03 

7.  4653E 

03 

5.  9722E 

03 

1.  4000E 

01 

6.  7703E 

03 

9.  6304E 

03 

7.  7033E 

03 

1.  4000E 

01 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

1.  OOOOE 

01 

6.  7703E 

03 

9.  6304E 

03 

7.  7033E 

03 

1.  OOOOE 

01 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

1 .  6000E 

01 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

1.  6000E 

01 

6.  7708E 

03 

9.  4184E 

03 

7.  5347E 

03 

1.  7000E 

01 

6.  7708E 

03 

9,  6304E 

03 

7.  7083E 

03 

1.  7000E 

01 

5.  2083E 

03 

7.  6S23E 

03 

6.  1453E 

03 

1.  8000 E 

01 

6.  7708E 

03 

9.  6304E 

03 

7.  70B3E 

03 

1.  OOOOE 

01 

6.  770SE 

03 

9.  6354E 

03 

7.  7083E 

03 

1 .  9000E 

01 

6.  7708E 

03 

9.  6334E 

03 

7.  7033E 

03 

1.  9000E 

01 

6.  7708E 

03 

9.  4184E 

03 

7.  5347E 

03 

2.  OOOOE 

01 

6.  770SE 

03 

9.  6304E 

03 

7.  70S3E 

03 

2.  OOOOE 

01 

6.  7703E 

03 

9.  6354E 

03 

7.  70S3E 

03 

2.  1000E 

01 

6.  7708E 

03 

9.  6354E 

03 

7.  7083E 

03 

2.  1000E 

01 

5.  2083E 

03 

7.  6823E 

03 

6.  1453E 

03 

2.  2000E 

01 

6.  770SE 

03 

9.  63S4E 

03 

7.  7083E 

03 

2.  2000E 

01 

6.  7708E 

03 

9.  4 1 S4E 

03 

7.  0347E 

03 

2.  3000E 

01 

6.  770SE 

03 

9.  6354E 

03 

7.  7083E 

03 

2.  3000E 

01 

6.  7703E 

03 

9.  6304E 

03 

7.  7033E 

03 

2.  4000E 

01 

6.  7708E 

03 

9.  6334E 

03 

7.  7083E 

03 

2.  4000E 

01 

6.  7708E 

03 

9.  6354E 

03 

7.  7083E 

03 

2.  5000E 

01 

6.  7708E 

03 

9.  6354E 

03 

7.  7083E 

03 

2.  5000E 

01 

5.  2083E 

03 

7.  4653E 

03 

5.  9722E 

03 

2.  6000E 

01 

6.  770BE 

03 

9.  6354E 

03 

7.  7033E 

03 

2.  6000E 

01 

6.  770SE 

03 

9.  6304E 

03 

7.  7033E 

03 

2.  7000E 

01 

6.  7708E 

03 

9.  6304E 

03 

7.  7083E 

03 

2.  700CE 

01 

6.  7708E 

03 

9.  6354E 

03 

7.  70S3E 

03 
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TIME  QO 


2. 

8000E 

01 

6.  7703E 

03 

2. 

8500E 

01 

6.  770SE 

03 

2. 

9000E 

01 

6.  7708E 

03 

2. 

9500E 

01 

5.  2083E 

03 

3. 

OOOOE 

01 

6.  770BE 

03 

3. 

OOOOE 

01 

6.  7708E 

03 

3. 

lOOOE 

01 

6.  7703E 

03 

3. 

1 500E 

01 

6.  7703E 

03 

3. 

2000E 

01 

6.  7708E 

03 

3. 

2500E 

01 

6.  7708E 

03 

3. 

30Q0E 

01 

6.  7708E 

03 

3. 

350CE 

01 

5.  2083E 

03 

3. 

4000E 

01 

6.  7708E 

03 

3. 

4500E 

01 

6.  7703E 

03 

3. 

5000E 

01 

6.  7708E 

03 

3. 

5000E 

01 

6.  770SE 

03 

3. 

6000E 

01 

6.  770SE 

03 

3. 

6500E 

01 

6.  7708E 

03 

3. 

7000E 

01 

6.  770SE 

03 

3. 

7500E 

Ol 

5.  2083E 

03 

3. 

8000E 

01 

6.  7708E 

03 

3. 

S300E 

01 

6.  7703E 

03 

3. 

9000E 

01 

6.  770SE 

03 

3. 

9500E 

01 

6.  770SE 

03 

4. 

OOOOE 

01 

6.  7708E 

03 

4. 

OOOOE 

01 

6.  7708E 

03 

4. 

lOOOE 

01 

6.  7708E 

03 

4. 

150OE 

01 

5.  2083E 

03 

4. 

2000E 

01 

6.  7708E 

03 

4. 

2500E 

01 

6.  770SE 

03 

4. 

3000E 

01 

6.  7708E 

03 

4. 

3500E 

01 

6.  7708E 

03 

4. 

4000E 

01 

6.  7708E 

03 

4. 

4500E 

01 

6.  7708E 

03 

4. 

OOOOE 

01 

6.  7708E 

03 

4. 

OOOOE 

01 

5.  2083E 

03 

4. 

600 OE 

01 

6.  770SE 

03 

4. 

6000E 

01 

6.  7708E 

03 

4. 

7000E 

01 

6.  7708E 

03 

4. 

7000E 

01 

6.  7708E 

03 

4. 

8000E 

01 

6.  770SE 

03 

QTOA 

QTO 

9.  6354E 

03 

7.  7083E 

03 

9.  4184E 

03 

7.  5347E 

03 

9.  6304E 

03 

7.  7083E 

03 

7.  6923E 

03 

6.  145SE 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6334E 

03 

7.  7083E 

03 

9.  4184E 

03 

7.  5347E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

7.  6823E 

03 

6.  1458E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  4184E 

03 

7.  5347E 

03 

9.  63S4E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6304E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7033E 

03 

7.  4633E 

03 

5.  9722E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6304E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  4184E 

03 

7.  5347E 

03 

9.  6304E 

03 

7.  7083E 

03 

7.  6823E 

03 

6.  145BE 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  4134E 

03 

7.  5347E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  6354E 

03 

7.  7083E 

03 

7.  6823E 

03 

6.  1458E 

03 

9.  6354E 

03 

7.  7083E 

03 

9.  4184E 

03 

7.  5347E 

03 

9.  6354E 

03 

7.  70B3E 

03 

9.  6304E 

03 

7.  70S3E 

03 

9.  6354E 

03 

7.  7083E 

03 
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VARIABLE 

MINIMUM 

T I  ME 

SNO 

2.  OOOOE 

01 

0.  OOOOE 

SNF 

3.  8100E 

00 

1.  9995E 

SCO 

2.  OOOOE 

02 

0.  OOOOE 

SCF 

1 .  6398E 

01 

1.  4000E 

S 

5.  OOOOE 

00 

8.  1562E 

X 

1 .  OOOOE 

02 

0.  OOOOE 

MAX  THUN  TIME! 


01 

7.  OOOOE 

01 

2.  OOOOE  00 

00 

5.  OOOOE 

01 

0. OOOOE -01 

01 

3.  OOOOE 

02 

2.  OOOOE  00 

01 

2.  OOOOE 

02 

0. OOOOE— 01 

00 

5.  OOOOE 

01 

0.  OOOOE— 01 

01 

1.  2369E 

02 

1.  OOOOE  01 
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INFLUENT  NITRATE  CONCENTRATION 


MINIMUM 
2.  5000E  01 


TIME 

SNO 

I 

0.  OOOOE- 

-01 

2.  5000E 

01 

4 

5.  OOOOE- 

-01 

2.  5000E 

01 

4 

1 .  OOOOE 

00 

2.  5000E 

01 

4 

1.  5000E 

00 

2.  5000E 

01 

4 

2.  OOOOE 

00 

7.  OOOOE 

01 

— 

2.  OOOOE 

00 

7.  OOOOE 

01 

— 

3.  OOOOE 

00 

7.  OOOOE 

01 

— 

3.  5000E 

00 

7.  OOOOE 

01 

— 

4.  OOOOE 

00 

7.  OOOOE 

01 

-- 

4.  OOOOE 

00 

2.  OOOOE 

01 

4 

0.  OOOOE 

00 

2.  OOOOE 

01 

4 

0.  OOOOE 

00 

2.  OOOOE 

01 

4 

6.  OOOOE 

00 

2.  OOOOE 

01 

4 

6.  OOOOE 

00 

2.  OOOOE 

01 

4 

7.  OOOOE 

00 

2.  OOOOE 

01 

4 

7.  OOOOE 

00 

2.  OOOOE 

01 

4 

8.  OOOOE 

00 

7.  OOOOE 

01 

— 

8.  OOOOE 

00  ' 

7.  OOOOE 

01 

-- 

9.  OOOOE 

00 

7.  OOOOE 

01 

-- 

9.  OOOOE 

00 

7.  OOOOE 

01 

— 

1 .  OOOOE 

01 

7.  OOOOE 

01 

— 

1.  OOOOE 

01 

2.  OOOOE 

01 

4 

1.  1000E 

01 

2.  OOOOE 

01 

4 

1.  1000E 

01 

2.  OOOOE 

01 

4 

1 .  2000E 

01 

2.  OOOOE 

01 

4 

1.  2000E 

01 

2.  OOOOE 

01 

4 

1.  3000E 

01 

2.  OOOOE 

01 

4 

1.  3000E 

01 

2.  OOOOE 

01 

•4 

1.  4000E 

01 

7.  OOOOE 

01 

— 

1.  4000E 

01 

7.  OOOOE 

01 

— 

1.  5000E 

01 

7.  OOOOE 

01 

— 

1.  OOOOE 

01 

7.  OOOOE 

01 

— 

1 .  6000E 

01 

7.  OOOOE 

01 

-■ 

1.  6000E 

01 

2.  OOOOE 

01 

4 

1.  7000E 

01 

2.  OOOOE 

01 

4 

1.  7000c 

01 

2.  OOOOE 

01 

4 

1 .  8000E 

01 

2.  OOOOE 

01 

4 

1 .  8000E 

01 

2.  OOOOE 

01 

4 

1.  900 OE 

01 

2.  OOOOE 

01 

4 

1.  9000E 

01 

2.  OOOOE 

01 

4 

2.  OOOOE 

01 

7.  OOOOE 

01 

- 

2.  OOOOE 

01 

7.  OOOOE 

01 

- 

2.  1000E 

01 

7.  OOOOE 

01 

- 

2.  1000E 

01 

7.  OOOOE 

01 

- 

2.  2000E 

01 

7.  OOOOE 

01 

- 

2.  2000E 

01 

2.  OOOOE 

01 

4 

2.  3000E 

01 

2.  OOOOE 

01 

4 

2.  3000E 

01 

2.  OOOOE 

01 

4 

2.  4000E 

01 

2.  OOOOE 

01 

4 

2.  4000E 

01 

2.  OOOOE 

01 

4 

2.  OOOOE 

01 

2.  OOOOE 

01 

4 

SNO  VERSUS  TIME  MAXIMUM 

7.  5000E  0 
I 


4 

■f 

4 

4 


4 

4 

4 

4 

4 


4 

4 

4 

4 

4 


4 

4 

4 

4 

4 
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INFLUENT  NITRATE  CONCENTRATION 


MINIMUM  SNO  VERSUS  TIME 

2.  5000E  01 

TIME  SNO  I 

2.  5500E  01  2.  SOOOE  01  + 

2.  6000E  01  7.  500 OE  01  - 

2.  6500E  01  7.  5000E  01  - 

2.  7000E  01  7.  5000E  01  - - - 

2.  7500E  01  7.  5000E  01  - 

2.  8000E  01  7.  5000E  01  - : - 

2.  B500E  01  2.  5000E  01  + 

2.  9000E  01  2.  5000E  01  + 

2.  9500E  01  2.  5000E  01  + 

3.  OOOOE  01  2.  5000E  01  + 

3.  0500E  01  2.  5000E  01  + 

3.  1000E  01  2.  5000E  01  + 

3.  1500E  01  2.  5000E  01  + 

3.  2000E  01  7.  5000E  01  - 

3.  2500E  01  7.  5000E  01  - 

3.  3000E  01  7.  5000E  01  - 

3.  3500E  01  7.  5000E  01  - 

3.  4000E  01  7.  5000E  01  - 

3.  4500E  01  2.  5000E  01  + 

3.  5000E  01  2.  5000E  01  + 

3.  5500E  01  2.  5000E  01  ■* 

3.  6000E  01  2.  5000E  01  + 

3.  6500E  01  2.  5000E  01  «• 

3.  7000E  01  2.  5000E  01  + 

3.  7500E  01  2.  5000E  01  + 

3.  8000E  01  7.  5000E  01  - 

3.  8500E  01  7.  SOOOE  01  - 

3.  9000E  01  7.  SOOOE  01  - 

3.  9500E  01  7.  SOOOE  01  - 

4.  OOOOE  01  7.  SOOOE  01  - 

4.  0500E  01  2.  SOOOE  01  + 

4.  1000E  01  2.  SOOOE  01  + 

4.  1500E  01  2.  SOOOE  01  + 

4.  2000E  01  2.  SOOOE  01  + 

4.  2500E  01  2.  SOOOE  01 

4.  3000E  01  2.  5000E  01  + 

4.  3500E  01  2.  SOOOE  01  + 

4.  4000E  01  7.  SOOOE  01  - 

4.  4500E  01  7.  SOOOE  01  - 

4.  SOOOE  01  7.  SOOOE  01  - 

4.  5500E  01  7.  SOOOE  01  - 

4.  6000E  01  7.  SOOOE  01  - 

4.  6500E  01  2.  SOOOE  01  + 

4.  7000E  01  2.  5000E  01  + 

4.  7500E  01  2.  SOOOE  01  + 

4.  8000E  01  2.  SOOOE  01  + 


MAXIMUM 
7.  5000E  01 
2 

- + 

- + 

- + 

- +  * 

- + 


+ 

+ 

■f 


- + 

- + 

- + 

- + 


EFFLUENT  NITRATE  CONCENTRATION 


TIME 

0.  0000E-01 
5.  0000E-01 
1.  OOOOE  00 

1.  OOOOE  00 

2.  OOOOE  00 

2.  OOOOE  00 

3.  OOOOE  00 

3.  OOOOE  00 

4.  OOOOE  00 

4.  OOOOE  00 

5.  OOOOE  00 

5.  OOOOE  00 

6.  OOOOE  00 

6.  OOOOE  00 

7.  OOOOE  00 

7.  OOOOE  00 

8.  OOOOE  00 

8.  OOOOE  00 

9.  OOOOE  00 
9.  OOOOE  00 
1.  OOOOE  01 
1.  OOOOE  01 
1.  1000E  01 
1.  1000E  01 
1.  2000E  01 
1.  2500E  01 
1.  3000E  01 
1.  3000E  01 
1.  4000E  01 
1.  4000E  01 
1.  OOOOE  01 
1.  OOOOE  01 
1.  6000E  01 
1.  6000E  01 
1.  7000E  01 
1.  7000E  01 
1.  8000E  01 
1.  8000E  01 
1.  9000E  01 

1.  9500E  01 

2.  OOOOE  01 
2.  OOOOE  01 
2.  1000E  01 
2.  1000E  01 
2.  2000E  01 
2.  2000E  01 
2.  3000E  01 
2.  3000E  01 
2.  4000E  01 
2.  4000E  01 
2.  OOOOE  01 


SNF 

0.  OOOOE 
9.  601  OE 
4.  0396E 
3.  S920E 

3.  8173E 
1.  4031 E 
1.  0492E 
1.  0701E 
1.  0731E 
6.  7026E 
0.  4609E 

4.  0107E 
3.  8322E 
0.  0432E 
0.  2164E 
0.  2412E 
0.  201 2E 
1.  4236E 
1.  0022E 
1.  1923E 
1.  14S3E 
6.  1413E 
0.  3734E 
0.  2636E 
0.  2479E 
0.  2407E 
0.  2404E 
3.  9834E 

3.  8289E 
1.  4033E 
1.  0492E 
1.  0701 E 
1.  0731 E 
6.  7026E 
0.  4609E 

4.  0107E 
3.  8322E 
0.  0432E 
0.  2164E 
0.  2412E 

5.  2012E 
1 .  4236E 
1.  0022E 
1.  1923E 
1.  14S3E 

6.  1413E 
0.  3734E 
5.  2636E 
0.  2479E 
5.  2407E 
0.  2404E 


MINIMUM 
3.  8100E  00 
I 


VERSUS  TIME 


MAXIMUM 
0.  OOOOE  01 
I 
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EFFLUENT  NITRATE  CONCENTRATION 


TIME 

SNF 

2.  5500E 

01 

3.  9834E 

2.  6000E 

01 

3.  8289E 

2.  6500E 

01 

1.  4033E 

2.  7000E 

01 

1.  5492E 

2.  7500E 

01 

1.  570 IE 

2.  8000E 

01 

1.  5731 E 

2.  8500E 

01 

6.  7526E 

2.  9000E 

01 

5.  4609E 

2.  9500E 

01 

4.  0107E 

3.  OOOOE 

01 

3.  8322E 

3.  0500E 

01 

5.  0432E 

3.  1000E 

01 

5.  2164E 

3.  1500E 

01 

5.  2412E 

3.  2000E 

01 

5.  2512E 

3.  2500E 

01 

1.  4236E 

3.  3000E 

01 

1.  5522E 

3.  3500E 

01 

1.  1923E 

3.  4000E 

01 

1.  1483E 

3.  4500E 

01 

6.  1413E 

3.  OOOOE 

01 

5.  3734E 

3.  OOOOE 

01 

5.  2636E 

3.  6000E 

01 

5.  2479E 

3.  6000E 

01 

5.  2457E 

3.  7000E 

01 

5.  2454E 

3.  7000E 

01 

3.  9834E 

3.  8000E 

01 

3.  828 9E 

3.  B500E 

01 

1 .  4033E 

3.  9000E 

01 

1.  5492E 

3.  9000E 

01 

1.  5701E 

4.  OOOOE 

01 

1.  573 IE 

4.  OOOOE 

01 

6.  7526E 

4.  1000E 

01 

5.  4609E 

4.  1500E 

01 

4.  0107E 

4.  2000E 

01 

3.  8322E 

4.  2000E 

01 

5.  0432E 

4.  3000E 

01 

5.  2164E 

4.  3000E 

01 

5.  2412E 

4.  4000E 

01 

5.  2512E 

4.  4500E 

01 

1.  4236E 

4.  OOOOE 

01 

1.  5522E 

4.  OOOOE 

01 

1.  1923E 

4.  600CE 

01 

1.  1483E 

4.  6000E 

01 

6.  141 3E 

4.  7000E 

01 

5.  3734E 

4.  7000E 

01 

5.  2636E 

4.  8000E 

01 

5.  2479E 

1-1  INI  HUH 
3.  8100E  00 


I 

00  + 

00  + 

0 1  - + 

0 1  - + 

01  - + 

01  - + 

00  - + 

00  -+ 

00  + 

00  + 

00  -+ 

00  -+ 

00  -+ 

00  -+ 

01  - + 

01  - + 

01  - + 

01  - + 

00  — + 

00  -+ 

00  -+ 

00  -  + 

00  -+ 

00  -+ 

00  + 

00  + 

01  - + 

01  - + 

01  - + 

01  - + 

00  - + 

00  ~+ 

00  + 

00  + 

00  -+ 

00  -+ 

00  -+ 

00  — 

01  - + 

CH  - + 

01  - + 

01  - + 

00  — + 

00  — 

00  -+ 

00  -+ 


BNF  VERSUS  TIME 


MAXIMUM 
5.  OOOOE  01 
I 


* 


s 
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T.  F.  INFLUENT  CONCENTRATION 


MINIMUM 
2.  OOOOE  02 


TIME 

SCO 

I 

0.  OOOOE - 

-01 

2.  5000E 

02 

+ 

5.  0QOOE- 

-01 

2.  OOOOE 

02 

t- 

1 .  OOOOE 

00 

2.  OOOOE 

02 

4 

1.  OOOOE 

00 

2.  OOOOE 

02 

f- 

2.  OOOOE 

00 

3.  OOOOE 

02 

— 

2.  OOOOE 

00 

3.  OOOOE 

02 

— 

3.  OOOOE 

00 

3.  OOOOE 

02 

— 

3.  OOOOE 

00 

3.  OOOOE 

02 

— 

4.  OOOOE 

00 

3.  OOOOE 

02 

4.  OOOOE 

00 

2.  OOOOE 

02 

f* 

0.  OOOOE 

00 

2.  OOOOE 

02 

•f 

0.  OOOOE 

00 

2.  OOOOE 

02 

4- 

6.  OOOOE 

00 

2.  OOOOE 

02 

4- 

6.  OOOOE 

00 

2.  OOOOE 

02 

4 

7.  OOOOE 

00 

2.  OOOOE 

02 

4- 

7.  OOOOE 

00 

2.  OOOOE 

02 

4- 

S.  OOOOE 

00 

3.  OOOOE 

02 

— 

S.  OOOOE 

00 

3.  OOOOE 

02 

-■ 

9.  OOOOE 

00 

3.  OOOOE 

02 

— 

9.  OOOOE 

00 

3.  OOOOE 

02 

— 

1 .  OOOOE 

01 

3.  OOOOE 

02 

— 

1.  OOOOE 

01 

2.  OOOOE 

02 

4- 

1.  1000E 

01 

2.  OOOOE 

02 

4- 

1.  1500E 

01 

2.  OOOOE 

02 

4- 

1.  2000E 

01 

2.  OOOOE 

02 

4- 

1.  250DE 

01 

2.  OOOOE 

02 

4- 

1.  3000E 

01 

2.  OOOOE 

02 

4- 

1.  3000E 

01 

2.  OOOOE 

02 

4- 

1.  4000E 

01 

3.  OOOOE 

02 

-• 

1.  4000E 

01 

3.  OOOOE 

02 

-■ 

1.  OOOOE 

01 

3.  OOOOE 

02 

— 

1.  OOOOE 

01 

3.  OOOOE 

02 

-■ 

1.  6000E 

01 

3.  OOOOE 

02 

— 

1.  6000E 

01 

2.  OOOOE 

02 

4- 

1.  7000E 

01 

2.  OOOOE 

02 

4 

1.  7000E 

01 

2.  OOOOE 

02 

4- 

1 .  8000E 

01 

2.  OOOOE 

02 

4- 

1.  BO0OE 

01 

2.  OOOOE 

02 

4- 

1.  9000E 

01 

2.  OOOOE 

02 

4- 

1 .  9500E 

01 

2.  OOOOE 

02 

4- 

2.  OOOOE 

01 

3.  OOOOE 

02 

-• 

2.  OOOOE 

01 

3.  OOOOE 

02 

-■ 

2.  1000E 

01 

3.  OOOOE 

02 

2.  1000E 

01 

3.  OOOOE 

02 

— 

2.  2000E 

01 

3.  OOOOE 

02 

— 

2.  2000E 

01 

2.  OOOOE 

02 

4- 

2.  3000E 

01 

2.  OOOOE 

02 

4- 

2.  3000E 

01 

2.  OOOOE 

02 

4- 

2.  4000E 

01 

2.  OOOOE 

02 

4* 

2.  4000E 

01 

2.  OOOOE 

02 

4- 

2.  OOOOE 

01 

2.  OOOOE 

02 

4- 

SCO  VERSUS  TIME  MAXIMUM 

3.  5000E  02 
I 


4- 

4- 

4- 


■4* 


+ 

4- 

4- 


4- 

4- 

4- 

4- 

4- 
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T.  F.  INFLUENT  CONCENTRATION 

MINIMUM 
2.  OOOOE  02 

TIME  SCO  I 


2. 

5500E 

01 

2.  OOOOE 

02 

4- 

2. 

6000E 

01 

3.  OOOOE 

02 

- 

2. 

6500E 

01 

3.  OOOOE 

02 

- 

2. 

7000E 

01 

3.  OOOOE 

02 

- 

2. 

75G0E 

01 

3  OOOOE 

02 

- 

2. 

8000E 

01 

3.  OOOOE 

02 

- 

2. 

B500E 

01 

2.  OOOOE 

02 

4- 

2. 

9000E 

01 

2.  OOOOE 

02 

4 

2. 

950QE 

01 

2.  OOOOE 

02 

4 

3. 

OOOOE 

01 

2.  OOOOE 

02 

4 

3. 

0500E 

01 

2.  OOOOE 

02 

4 

3. 

1 OOOE 

01 

2.  OOOOE 

02 

4 

3. 

1500E 

01 

2.  OOOOE 

02 

4 

3. 

2000E 

01 

3.  OOOOE 

02 

- 

3. 

2500E 

01 

3.  OOOOE 

02 

- 

3. 

3000E 

01 

3.  OOOOE 

02 

- 

3. 

3500E 

01 

3.  OOOOE 

02 

- 

3. 

4000E 

01 

3.  OOOOE 

02 

- 

3. 

4500E 

01 

2.  OOOOE 

02 

4- 

3. 

5000E 

01 

2.  OOOOE 

02 

4 

3. 

5500E 

01 

2.  OOOOE 

02 

4- 

3. 

6000c 

01 

2.  OOOOE 

02 

4 

3. 

6500E 

01 

2.  OOOOE 

02 

4 

3. 

7000E 

01 

2.  OOOOE 

02 

4 

3. 

7500E 

01 

2.  OOOOE 

02 

4- 

3. 

BOOOE 

01 

3.  OOOOE 

02 

- 

3. 

S500E 

01 

3.  OOOOE 

02 

- 

3. 

9000E 

01 

3.  OOOOE 

02 

- 

3. 

9500E 

01 

3.  OOOOE 

02 

- 

4. 

OOOOE 

01 

3.  OOOOE 

02 

- 

4. 

050GE 

01 

2.  OOOOE 

02 

4- 

4. 

1000E 

01 

2.  OOOOE 

02 

4- 

4. 

1 500E 

01 

2.  OOOOE 

02 

-4- 

4. 

2000E 

01 

2.  OOOOE 

02 

4- 

4. 

2500E 

01 

2.  OOOOE 

02 

4- 

4. 

3000E 

01 

2.  OOOOE 

02 

•4- 

4. 

3500E 

01 

2.  OOOOE 

02 

4- 

4. 

4000E 

01 

3.  OOOOE 

02 

- 

4. 

4500E 

01 

3.  OOOOE 

02 

- 

4. 

5000E 

01 

3.  OOOOE 

02 

- 

4. 

5  OOOE 

01 

3.  OOOOE 

02 

- 

4. 

6000E 

01 

3-  OOOOE 

02 

- 

4. 

6  OOOE 

01 

2.  OOOOE 

02 

4- 

4. 

7000E 

01 

2.  OOOOE 

02 

4- 

4. 

7  OOOE 

01 

2.  OOOOE 

02 

4- 

4. 

BOOOE 

01 

2.  OOOOE 

02 

4- 

SCO  VERSUS  TIME  MAXIMUM 

3.  5000E  02 
1 

- + 

- + 

- + 

- +  , 

- + 
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T.  F.  EFFLUENT  CONCENTRATION 


MINIMUM  SCF  VERSUS  TIME 


1. 

6398E  i 

TIME 

SCF 

i 

0.  OOOOE- 

-01 

2.  OOOOE 

02 

— 

5.  OOOOE- 

-01 

7.  1  54 OE 

01 

— 

1 .  OOOOE 

00 

2.  9192E 

01 

- 4- 

1.  OOOOE 

00 

1.  9140E 

01 

4 

2.  OOOOE 

00 

1.  676 IE 

01 

4* 

2.  OOOOE 

00 

2.  7803E 

01 

— 4. 

3.  OOOOE 

00 

3.  080 IE 

01 

- 4. 

3.  OOOOE 

00 

3.  161  OE 

01 

- 4 

4.  OOOOE 

00 

3.  1836E 

01 

- + 

4.  OOOOE 

00 

2.  4724E 

01 

— 4 

0.  OOOOE 

00 

2.  2842E 

01 

—4* 

5.  OOOOE 

00 

1.  8169E 

01 

4 

6.  OOOOE 

00 

1.  7001 E 

01 

4 

6.  OOOOE 

00 

2.  123SE 

01 

— 4 

7.  OOOOE 

00 

2.  2375E 

01 

—4 

7.  OOOOE 

00 

2.  2196E 

01 

-•4 

B.  OOOOE 

00 

2.  2130E 

01 

—4 

B.  OOOOE 

00 

2.  926 IE 

01 

- 4 

9.  OOOOE 

00 

3.  1 1 97E 

01 

- 4 

9.  OOOOE 

00 

2.  024SE 

01 

—  4 

1.  OOOOE 

01 

2.  3S26E 

01 

— 4 

1.  OOOOE 

01 

2.  2569E 

01 

—  4 

1.  1000E 

01 

2.  2263E 

01 

-4 

1.  1000E 

01 

2.  2603E 

01 

— 4 

1.  2000E 

01 

2.  2759E 

01 

— 4 

1.  2000E 

01 

2.  278BE 

01 

— 4 

1.  300 OE 

01 

2.  2796E 

01 

—  4 

1.  3500E 

01 

1.  7634E 

01 

4 

1 . 4000E 

01 

1 .  6404E 

01 

4- 

1.  4000E 

01 

2.  7706E 

01 

- 4 

1.  OOOOE 

01 

3.  0774E 

01 

- 4 

1.  OOOOE 

01 

3.  1600E 

01 

- 4 

I.  6000E 

01 

3.  1834E 

01 

- 4 

1.  6500E 

01 

2.  4723E 

01 

— 4 

1.  7000E 

01 

2.  2842E 

01 

-4 

1.  7000E 

01 

1.  8169E 

01 

4 

1.  8000E 

01 

1.  7051E 

01 

4 

1.  8000E 

01 

2.  1 238E 

01 

-4 

1.  9000E 

01 

2.  237 OE 

01 

-4 

1.  9000E 

01 

2.  2196E 

01 

-4 

2.  OOOOE 

01 

2.  2130E 

01 

—  4 

2.  OOOOE 

01 

2.  9261 E 

01 

- 4 

2.  1000E 

01 

3.  1197E 

01 

- H- 

2.  1000E 

01 

2.  5248E 

01 

•—4 

2.  2000E 

01 

2.  3S26E 

01 

— 4 

2.  2000E 

01 

2.  2069E 

01 

—  4 

2.  3000E 

01 

2.  2263E 

01 

—  4 

2.  3000E 

01 

2.  2653E 

01 

-4 

2.  4000E 

01 

2.  2759E 

01 

~4 

2.  4000E 

01 

2.  2788E 

01 

-4 

2.  OOOOE 

01 

2.  2796E 

01 

—4 

MAXIMUM 
2.  OOOOE  02 
I 
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T  F.  EFFLUENT  CONCENTRATION 


MINIMUM 

SCF 

VERSUS  TIME 

MAX  I HUM 

1.6398E  01 

p.  5oooe  o; 

TIME 

SCF 

I 

I 

2. 

55Q0E 

01 

1.  7634  E 

01 

4 

2. 

6Q00E 

01 

1.  6404E 

01 

4- 

2. 

650DE 

01 

2.  7706E 

01 

- 1- 

2. 

7000E 

01 

3.  0774F 

01 

- 4- 

2. 

7500E 

01 

3.  1608E 

01 

- 4- 

2. 

8000E 

01 

3.  1834E 

01 

- 4- 

2. 

8500E 

01 

2.  4723E 

01 

— • 4 

2. 

9000E 

01 

2.  2842E 

01 

-  4- 

2. 

9500E 

01 

1.  S169E 

01 

4 

3. 

OOOOE 

01 

1.  7051 E 

01 

4- 

3. 

0500E 

01 

2.  1238E 

01 

~4 

3. 

1 OOOE 

01 

2.  2375E 

01 

— 4 

3. 

1 50  OE 

01 

2.  2196E 

01 

-4- 

3. 

2000E 

01 

2.  2130E 

01 

—4- 

3. 

2500E 

01 

2.  926 IE 

01 

- + 

3. 

3000E 

01 

3.  1197E 

01 

- 4 

3. 

3500E 

01 

2.  S248E 

01 

— 4 

3. 

4000E 

01 

2.  3S26E 

01 

-4- 

3. 

4500E 

01 

2.  2569E 

01 

— f- 

3. 

500 OE 

01 

2.  2263E 

01 

-4- 

3. 

5500E 

01 

2.  2653E 

01 

— f 

3. 

6000E 

01 

2.  2759E 

01 

-  4- 

3. 

6500E 

01 

2.  2788E 

01 

-4 

3. 

7000E 

01 

2.  2796E 

01 

-4- 

3. 

7500E 

01 

1.  7634E 

01 

4 

3. 

8000E 

01 

1.  640 'IE 

01 

4 

3. 

8500E 

01 

2.  7706E 

01 

- 4 

3. 

9000E 

01 

3.  0774E 

01 

- 4 

3. 

9500E 

01 

3.  1608E 

01 

- 4 

4. 

OOOOE 

01 

3.  1S34E 

01 

- 4- 

4.  0500E 

01 

2.  4723E 

01 

— 4- 

4. 

1  OOOE 

01 

2.  2842E 

01 

-4 

4. 

1 50CE 

01 

1.  8169E 

01 

4* 

4. 

2000E 

01 

1 . 7 0  5 1 E 

01 

4 

4. 

2500E 

01 

2.  1238E 

01 

-4- 

4. 

3000E 

01 

2.  2375E 

01 

-4 

4. 

3500E 

01 

2.  2196E 

01 

—4- 

4. 

4000E 

01 

2.  2130E 

01 

— f- 

4. 

4500E 

01 

2.  926 IE 

01 

- 4- 

4. 

5000E 

01 

3.  I197E 

01 

- 4- 

4. 

5500E 

01 

2.  5248E 

01 

—  4- 

4. 

6OO0E 

01 

2.  3826E 

01 

-4- 

4. 

6500E 

01 

2.  2569E 

01 

— (- 

4. 

7000E 

01 

2.  226 3E 

01 

—  4" 

4. 

7500E 

01 

2.  2653E 

01 

— 4 

4. 

8000E 

01 

2.  2759E 

01 

-4- 

A.  S.  EFFLUENT  CONCENTRATION 


TINE 

0.  OOOOE- 

nnn.'ic. 

-01 

s 

5.  OOOOE 

1' 

5. 

01 
n  i 

IINIMUH  S  VERSUS  TIME  MAXIMUM 

OOOOE  00  0.  OOOOE  0] 

I  1 

1.  0000E 

00 

4.  1637E 

01 

1.  OOOOE 

00 

3.  7 1  2 1 F 

0.1 

- - - +■ 

2.  OOOOE 

00 

3.  2770E 

01 

2.  OOOOE 

00 

2.  S755E 

01 

- - - + 

3.  OOOOE 

00 

2.  5145E 

01 

- + 

3.  OOOOE 

00 

2.  1SS5E 

01 

- +• 

4.  OOOOE 

00 

1.  896 IE 

01 

- + 

4.  OOOOE 

00 

1.  6305E 

01 

- 4. 

0.  OOOOE 

00 

1 .  3923E 

01 

- + 

0.  OOOOE 

00 

1.  1793E 

01 

- + 

6.  OOOOE 

00 

9.  934 IE 

00 

- 4. 

6.  OOOOE 

00 

8.  397 IE 

00 

— - f" 

7.  OOOOE 

00 

7.  1307E 

00 

- 4. 

7.  OOOOE 

00 

6.  0725E 

00 

-•4- 

8.  OOOOE 

00 

5.  1919E 

00 

i 

8.  OOOOE 

00 

5.  OOOOE 

00 

4* 

9.  OOOOE 

00 

5.  OOOOE 

00 

4 

9.  OOOOE 

00 

5.  OOOOE 

00 

4 

1.  OOOOE 

01 

5.  OOOOE 

00 

4 

1.  OOOOE 

01 

5.  OOOOE 

00 

4- 

1.  1000E 

01 

5.  OOOOE 

00 

4 

1.  1500E 

01 

5.  OOOOE 

00 

4' 

1.  2000E 

01 

5.  OOOOE 

00 

4 

1.  2000E 

01 

5.  OOOOE 

00 

4* 

1.  3000E 

01 

5.  OOOOE 

00 

4- 

1.  3500E 

01 

5.  OOOOE 

00 

4 

1.  4000E 

01 

5.  OOOOE 

00 

4- 

1.  4000E 

01 

0.  OOOOE 

00 

4- 

1.  OOOOE 

01 

5.  OOOOE 

00 

4- 

1.  OOOOE 

01 

5.  OOOOE 

00 

4- 

1.  6000E 

01 

5.  OOOOE 

00 

4- 

1.  6500E 

01 

5.  OOOOE 

00 

4- 

1.  7000E 

01 

5.  OOOOE 

00 

4* 

1.  7000E 

01 

5.  OOOOE 

00 

4* 

1.  8000E 

01 

5.  OOOOE 

00 

4. 

1 . SOOOE 

01 

5.  OOOOE 

00 

4* 

1 . 9000E 

01 

5.  OOOOE 

00 

4- 

1.  9000E 

01 

5.  OOOOE 

00 

4- 

2.  OOOOE 

01 

5.  OOOOE 

00 

4* 

2.  OOOOE 

01 

5.  OOOOE 

00 

4- 

2.  1000E 

01 

5.  OOOOE 

00 

4- 

2.  1000E 

01 

5.  OOOOE 

00 

4- 

2.  2000E 

01 

5.  OOOOE 

00 

4- 

2.  2000E 

01 

5.  OOOOE 

00 

4* 

2.  3000E 

01 

5.  OOOOE 

00 

4* 

2.  SOOOE 

01 

5.  OOOOE 

00 

4- 

2.  4000E 

01 

5.  OOOOE 

00 

4- 

2.  4000E 

01 

0.  OOOOE 

00 

4* 

2.  OOOOE 

01 

5.  OOOOE 

00 

4- 

A.  S.  EFFLUENT  CONCENTRATION 


MINIMUM 

S  VERSUS  TIME 

MAXIMUM 

0.  OOOOE  00 

5.  OOOOE  01 

TIME 

S 

I 

I 

2. 

OOOOE 

01 

5.  OOOOE 

00 

4 

2. 

6000E 

01 

5.  OOOOE 

00 

4* 

2. 

6500E 

01 

5.  OOOOE 

00 

4 

2. 

7000E 

01 

5.  OOOOE 

00 

4 

2. 

7500E 

01 

5.  OOOOE 

00 

4 

2. 

BOOOE 

01 

5.  OOOOE 

00 

4 

2. 

8500E 

01 

5.  OOOOE 

00 

4 

2. 

9000E 

01 

5.  OOOOE 

00 

4- 

2. 

9500E 

01 

5.  OOOOE 

00 

4* 

3. 

OOOOE 

01 

5.  OOOOE 

00 

4 

3. 

OOOOE 

01 

5.  OOOOE 

00 

4- 

3. 

1000E 

01 

5.  OOOOE 

00 

4- 

3. 

1 50  OE 

01 

5.  OOOOE 

00 

4* 

3. 

2000E 

01 

5.  OOOOE 

00 

4 

3. 

2500E 

01 

5.  OOOOE 

00 

4 

3. 

3000E 

01 

5.  OOOOE 

00 

4 

3. 

3500E 

01 

5.  OOOOE 

00 

4 

3. 

4000E 

01 

5.  OOOOE 

00 

4* 

3. 

4500E 

01 

5.  OOOOE 

00 

4- 

3. 

OOOOE 

01 

5.  OOOOE 

00 

4* 

3. 

OOOOE 

01 

5.  OOOOE 

00 

4* 

3. 

6000E 

01 

5.  OOOOE 

00 

4* 

3. 

6500E 

01 

5.  OOOOE 

00 

4 

3. 

7000 E 

01 

5.  OOOOE 

00 

4- 

3. 

7000E 

01 

5.  OOOOE 

00 

4- 

3. 

OOOOE 

01 

5.  OOOOE 

00 

4 

3. 

OOOOE 

01 

5.  OOOOE 

00 

4 

3. 

9000E 

01 

5.  OOOOE 

00 

4 

3. 

9500E 

01 

5.  OOOOE 

00 

4 

4. 

OOOOE 

01 

5.  OOOOE 

00 

4 

4. 

OOOOE 

01 

5.  OOOOE 

00 

4- 

4. 

1 OOOE 

01 

5.  OOOOE 

00 

4- 

4. 

1  OOOE 

01 

5.  OOOOE 

00 

4 

4. 

2000E 

01 

5.  OOOOE 

00 

4 

4. 

2500E 

01 

5.  OOOOE 

00 

4- 

4. 

3000E 

01 

5.  OOOOE 

00 

4 

4. 

3  OOOE 

01 

5.  OOOOE 

00 

4* 

4. 

4  OOOE 

01 

5.  OOOOE 

00 

•4 

4. 

4  OOOE 

01 

5.  OOOOE 

00 

4 

4. 

OOOOE 

01 

5.  OOOOE 

00 

4- 

4. 

OOOOE 

01 

5.  OOOOE 

00 

4 

4. 

6000E 

01 

5.  OOOOE 

00 

4- 

4. 

6500E 

01 

5.  OOOOE 

00 

4* 

4. 

7  OOOE 

01 

5.  OOOOE 

00 

4 

4. 

7  OOOE 

01 

5.  OOOOE 

00 

4- 

4. 

BOOOE 

01 

5.  OOOOE 

00 

4- 
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ORGAN  I BM  COIMCENTR  AT  I  ON 


TIME 


MINIMUM 
1.  OOOOE  02 


O.  OOOOE— 01 

5.  OOOOE -01 

1.  OOOOE  00 
1  5000E  00 

2.  OOOOE  00 

2.  OOOOE  00 

3.  OOOOE  00 

3.  OOOOE  00 

4.  OOOOE  00 

4.  OOOOE  00 
0.  OOOOE  00 

5.  OOOOE  00 

6.  OOOOE  00 

6.  OOOOE  00 

7.  OOOOE  OO 

7.  OOOOE  00 

8.  OOOOE  00 
3.  OOOOE  00 

9.  OOOOE  00 
9.  OOOOE  00 
1.  OOOOE  01 
1.  OOOOE  01 
1.  1000E  01 
1.  1500E  01 
1.  2000E  01 
1.  2000E  01 
1.  3000E  01 
1.  3500E  01 
1.  4000E  01 
1.  45  ODE  01 
1.  OOOOE  01 
1.  5500E  01 
1.  6000E  01 
1.  6000E  01 
1.  7000E  01 
1.  7500E  01 
1.  8000E  01 
1.  S500E  01 
1.  9000E  01 

1.  9500E  01 

2.  OOOOE  01 
2.  0500E  01 
2.  1000E  01 
2.  1500E  01 
2.  2000E  01 
2.  2500E  01 
2.  3000E  01 
2.  3500E  01 
2.  4000E  01 
2.  4500E  01 
2.  5000E  01 


1 .  OOOOE  02 
1 .  026BE  02 
1 .  0527E  02 
1  0772E  02 
1.  1001 E  02 
1.  1210E  02 
1.  140 IE  02 
1.  1572E  02 
1.  1725E  02 
1.  1B53E  02 
1.  1973E  02 
1 .  2073E  02 
1.  2155E  02 
1.22.17E  02 
1.  226 SE  02 
1.  230 IE  02 
1.  2327E  02 
1.  234 4E  02 
1.  2355E  02 
1.  236 5E  02 
1 .  2369E  02 
1 . 2364E  02 
1.  2356E  02 
1.  2344E  02 
1.  2331E  02 
1.  2315E  02 
1 .  2298E  02 
1 .  2236E  02 
1 .  227 IF.  02 
1.  2250E  02 
1 .  2229E  02 
1 .  2209E  02 
1.  2191E  02 
1.  2173E  02 
1.  2156E  02 
1.  2 1 4 1 E  02 
1.  2125E  02 
1.2104E  02 
1 .  2033E  02 
1 .  2063E  02 
1.  204  BE  02 
1.  2022E  02 
1 .  2003E  02 
1.  1999E  02 
1.  1975E  02 
1.  1957E  02 
1.  1938E  02 
1.  1920E  02 
1.  190 IE  02 
1.  1S33E  02 
1.  1865E  02 


X  VERSUS  TIME  MAXIMUM 

1 .  2369E  02 
I 
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ru  ru  ru  ru  ru 


ORGAN I BM  COWCCNTR AT  I  ON 


TINE 

.  5500E  Ot 
.  600 or;  o  1 
.  6500E  01 
.  7000c  0  I 
.  7500E  01 
2.  SOOOE  01 
2.  8500E  01 
2.  9000E  01 

2.  9500E  01 

3.  OOOOE  01 
3.  OOOOE  01 
3.  1000E  01 
3.  1500E  01 
3.  2000E  01 
3.  2D00E  01 
3.  3000E  01 
3.  3500E  01 
3.  4000E  01 
3.  4500E  01 
3.  SOOOE  01 
3.  5500E  01 
3.  6000E  01 
3.  6500E  01 
3.  7000E  01 
3.  7500E  01 
3.  SOOOE  01 
3.  B500E  01 
3.  9000E  01 

3.  9500E  01 

4.  OOOOE  01 
4.  0500E  01 
4. 1000E  01 
4.  1500E  01 
4.  2000E  01 
4.  2500E  01 
4.  3000E  01 
4.  350CE  01 
4. 4000E  01 
4.  4500E  01 
4.  SOOOE  01 
4.  5500E  01 
4.  6000E  01 
4.  6500E  01 
4.  7000E  01 
4.  7500E  01 
4.  8000E  01 


MINIMUM  X  VERSUS  TIME  MAXIMUM 

1.  OOOOE  02  1.  2369E  02 

X  I  I 

1.  1QS2E  02  - + 

1.  1337E  02  - + 

1.  lSIV'E  02  - + 

1.  17S2E  02  +  1 

1.  1767E  02  + 

1.  17S3E  02  + 

1.  1739F  02  + 

1  1720E  02  + 

1.  1716E  02  + 

1.  1699E  02  + 

1.  1682E  02  + 

1.  1666E  02  + 

1.  .16S0E  02  + 

1.  1634E  02  + 

1.  1620E  02  + 

1.  1 60 9 E  02  + 

1.  1S9SE  02  + 

1.  1S84E  02  + 

1.  1S70E  02  + 

1.  15S6E  02  + 

1.  1542E  02  + 

1.  1  52 BE  02  + 

1.  1S14E  02  + 

1.  1S04E  02  + 

1.  1492E  02  + 

1.  1477E  02  + 

1.  1463E  02  + 

] .  1  4S0E  02  + 

1.  1439E  02  + 

1.  1429E  02  + 

1.  141 SE  02  - + 

1.  141 OE  02  + 

1.  140 IE  02  + 

1.  13SSE  02  + 

1.  1  37 SE  02  + 

1.  1362E  02  + 

1.  13S0E  02  + 

1.  1338E  02  + 

1.  1327E  02  + 

1.  1319E  02  + 

1.  131  IE  02  + 

1.  1300E  02  +  * 

1.  1 29 OE  02  + 

1.  1279E  02  - + 

1.  1268E  02  - + 
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APPENDIX  B 


VARIABLE  LIST 
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The  program  written  in  BASIC  uses  the  same  model  as  the  CSMP  program.  The 
programs  use  Runge-Kutta  as  a  method  of  integration. 

The  variables  used  in  the  two  programs  are  described  in  the  following 


pages;  the  BASIC  program  is  arranged  as 
Line  # 

100  -  900 

1000  -  1999 

2000  -  2999 

3000  -  3999 


follows: 

Function 

constants,  initial  conditions, 
pulses  set 

anoxic  tower  with  lines 
1500-1599  as  subroutine 

trickling  filter  with  lines 
2500-2999  as  subroutine 

activated  sludge  with  lines 
3500-3999  as  subroutine 


VARIABLE  LIST 


ANOXIC  TOWER 


CSMP 

BASIC 

A 

Cross  sectional  area  (ft) 

A 

ZLAM 

Rate  constant  (gal/min-f t^) 

ZLAM 

U 

Packing  coefficient 

U 

R 

Effectiveness  factor 

R 

Z 

Column  depth  (ft) 

Z 

VOL 

Volume  (ft^) 

VOL 

QINZ 

Steady-state  flow  (mgd) 

QZ 

QIN1 

Steady-state  flow  (gpm) 

QA 

QIN3 

Additional  peak  flow  (mgd) 

Q3 

QINF 

Additional  peak  flow  (gpm) 

QB 

QO 

Flow  thru  anoxic  tower  (gpm) 

QO 

QZ 

Conversion  of  QO  to  ft^/hr 

SNI 

Initial  condition  -  nitrates  (mg/L) 

SNI 

SNINF 

Influent  nitrates  (mg/L) 

SN 

SNO 

Influent  nitrates  (pulse) 

NS 

SNI 

Peak  flow  -  nitrates  (mg/L) 

SI 

SNF 

Effluent  -  nitrate  (mg/L) 

FS 

TRICKLING  FILTER 

CA 

Cross  sectional  area  (ft) 

CA 

CLAM 

Rate  constant  (gal/min  -  ft^) 

CLAM 

ZT 

Depth  (ft) 

ZT 

W 

Packing  coefficient 

W 

B 

Effectiveness  factor 

B 

CVOL 

Volume  (ft^) 

CVOL 

QIC2 

Constant  flow  rate  (mgd) 

Q4 

QIC1 

Constant  flow  (gpm) 

QC 

QIC3 

Pulse  flow  (mgd) 

Q5 

QICF 

Pulse  flow  (gpm) 

QD 

QTO 

Flow  thru  trickling  filter  (gpm) 

QT 

Ql 

Conversion  of  QTO  to  ft^/hr 

SCI 

Initial  condition  -  carbon  (mg/L) 

Cl 

SCI 

Pulse  cone.  -  carbon  (mg/L) 

Cl 

SCINF 

Carbonaceous  substrate  influent  (mg/L) 

CF 

SCO 

Pulse-carbonaceous  substrate  influent 

SO 
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VARIABLE  LIST 


ACTIVATED  SLUDGE 


CSMP 

BASIC 

K 

Rate  constant  ( 1/day) 

K 

K1 

K  ( 1/hr) 

KS 

Half-velocity  constant  (mg/L) 

KS 

KD 

Endogenous  decay  coefficient  (1/day) 

KD 

K2 

KD  (1/hr) 

YG 

Yield  coefficient 

YG 

AV0L1 

Volume  (ft^) 

AVOL 

QR 

Set  recirculation 

QR 

XI 

Initial  condition  on  organism  (mg/L) 

XI 

SC 

Initial  condition  on  substrate  (mg/L) 

SC 

XC 

Startup  condition  on  organism  (mg/L) 

XC 

XO 

Organism  concentration  with  recycle  (mg/L) 

XO 

X 

Organism  concentration  (mg/L) 

X 

S 

Substrate  concentration  (mg/L) 

S 

QTOA 

Influent  flow  to  activated  sludge  (gpm) 

QE 
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ADDITIONAL  VARIABLES  -  BASIC: 


M8  =  Counter  (outputs  every  5th  valve) 

T,  T2,  T3  =  Time  for  anoxic  tower,  trickling  filter,  activated  sludge, 
respectively 
Kl,  K2 ,  K3 ,  K4 

K5 ,  K6,  K7 ,  K8  Valves  for  Runge-Kutta 
S8 ,  X8 

T4,  T5,  T6  =  Counter  for  time  array 

EQUATIONS  DESCRIBING  PROCESS: 

ANOXIC  TOWER: 

DN  =  U1  -  [EXP( U2) ]  *  FS 
where 

U1  =  (QO/VOL)  *  (NS-FS) 

UZ  =  [-ZLAM  *  (Z*R)]/[(QO/A)*U] 

TRICKLING  FILTER: 

DC  =  U3  -  [  EXP(  U4)  ]  *  SF 
where 

U3  =  (QT/CVOL)  *  (SO  -  SF) 

U4  =  [-CLAM*  [Z'TB)]/[(QT/CA)*W] 

ACTIVATED  SLUDGE: 

X(I)  =  (US  -  U6)  *  H 
where 

U5  =  (QE/AVOL)  *  (X0  -  X)  *  1440 

U6  =  [((K  *  YG)  *  (S/CS  +  KS))  -  KD  *  X)] 

H  =  Increment  for  Runge-Kutta 

S(  I)  -  H  *  [((QE/AVOL)  *  (SF  ( T3)  -  S))  *  1440  -  (K  *  (S/CS  +  KS))  *  X)  ] 


ADDITIONAL  VARIABLES  NEEDED  FOR  BASIC  PROGRAM: 


PI: 

Time  where  pulse  begins  (for 

TRIG1  in  CSMP) 

P2 : 

Time  spacing  between  pulses 

(for  TRIG1) 

P3 : 

Length  of  pulse  (for  TRIG1) 

P4 : 

Time  where  pulse  beings  (for 

TRIG2) 

P5 : 

Time  spacing  between  pulses 

P6 : 

Length  of  pulse 

P7: 

Begin  pulse  (TRIG  3) 

P8 : 

Spacing 

P9 : 

Length 

P0 : 

Begin  pulse  (TRIG4) 

PI: 

Spacing 

PJ: 

Length 

G,G1:  To  keep  pulses  moving  with  time: 


A 


o 

CO 

r> 

C/3 


Width 


r — 


Start  i 

i- - 4- 


Spacing 


G,G1=0 


- •*  ;  4  G,G1  =  1 


Time  — + 


APPENDIX  C 

BASIC  MODEL 
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20  LPR I  NT : LPR I N  T : LPR I NT 

30  LPR I NT  CHR$ (14)"  INDUSTRIAL  WASTE  TREATMENT  PLANT" 

40  LPR I NT: LPR I NT 

50  LPR I NT  CHR$(14)"  HOLS TON  AAP" 

60  LPR I NT : LPR I NT : LPR I NT 
65  LPR  I  NT LPR  I  NT ::  LPR  I  NT 

70  LPR I NT  TAB (15); "HOW  MANY  HOURS  WOULD  YOU  LIKE  TO  SEE  ?" 

72  PR I NT "HOW  MANY  HOURS  WOULD  YOU  LIKE  TO  SEE?" 

75  INPUT  "CHOOSE  UP  TO  48";LP 
77  LPR  I  NT  TAB  (15)  ;  L..P 
SO  LPR I NT ; LPR I NT : LPR I NT 

82  INPUT  "INPUT  FLOW  FROM  AREA  B  (MOD)  " ; QB 

83  LPR I NT: LPR I NT 

84  LPR I NT  TAB  015) 5 "FLOW  FROM  AREA  B  (MOD)  IS 
"QB: LPR I  NT: LPR I NT 

85  INPUT  "INPUT  FLOW  FROM  AREA  A  (MOD)  " ? QD 

86  LPR I NT  TAB (15); "FLOW  FROM  AREA  A  (MOD)  IS 
"QD: LPR I NT: LPR I NT 

87  LPR I NT  TAB < 1 5) 5 " ENTER  VALUE  FOR  COD  ( mg / L )  I N  TR I CKL I NG 
F I LTER  I NF  L DENT " : I NPUT  CF : LPR I NT  TAB ( 1 5 ) 5  CF 

88  LPR I NT  TAB (15) 5 "DO  YOU  WANT  THE  TRICKLING  FILTER  EFFLUENT 
TO  BE  THE": LPR I NT  TAB (15) 5 "ACTIVATED  SLUDGE  INFLUENT  ?  (Y  or 
N)": INPUT  X* 

89  LPR I NT  TAB ( 15) 5  X$: LPR I NT: LPR I NT 

90  D I M  GO  ( 500 )  ,  QT  <  500 )  ,  S  ( 500 )  ,,  X  ( 500 )  ,,  SF  ( 500 ) 

98  DIM  T (500) 

.100  REM  ...CONSTANTS 
110  REM  „ „ .ANOXIC  TOWERS 
120  A  6655!  :  REM  ..  .AREA  (C.S.  ) 

130  ZLAM  *'  1.45: REM  . . .  RATE  CONSTANT  (CF/CF-HR) 

135  IJ  *  .5:  REM  ...PACKING  COEFF. 

140  R  =  .6:  REM  . EFFECTIVENESS  FACTOR  FOR  TOWER 

150  Z  =  10!:  REM  ...DEPTH  OF  COLUMN  MEDIA 
160  VOL.  *  66550 1  :  REM  ...VOLUME 

170  Q2  =  15!:  REM  ...CONSTANT  FLOW  RATE 

171  REM  .  .  .  FLOW  RATES  ARE  MOD 

180  Q3  -  4.5  : REM  . PEAK  NITRATES 

190  SN  —  25 ! : REM  ...  PEAK  NI i RA I ES 
200  SI  ~  10!: REM  ...INFLUENT  NITRATES 
210  SNI  ••=  25!:  REM  .  .  .  INIT.  COND.  -NI  T RATES 
220  REM  ...TRICKLING  FILTER 
230  CA  ~  10000: REM  ...AREA  (C.S.) 

240  CLAM  =  .33:  REM  ...RATE  CONSTANT  FOR  TRICKLING  FILTER 

250  ZT  =  26  !  :  REM  .  .  DEPTH 

260  W  =  .5;  REM  ...PACKING  COEFF. 

270  B  -  .67: REM  ...EFFECTIVENESS  FACTOR 
280  CVOL  =  260000! 

290  Q4  «  1.1:  REM  ...CONS.  FLOW  ( MGD ) 

300  Q5  =  .25:  REM  ...ADDITIONAL  PEAK  FLOW 

310  REM  ...  CF  ...  SUBSTRATE  C  I NFL . " 

320  Cl  =  100 

330  Cl  =  50  : REM  ...INIT.  COND.  (MG/L) 

340  REM  . . .ACTIVATED  SLUDGE 

350  K  -  5: REM  ...  RATE  CONSTANT  ( 1 / DAY ) 

351  K  ~  K/24:  REM  ...CONVERT  TO  (1/HR.) 
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360 

370 

371 
380 
390 
4 1 0 
420 
430 
440 
450 
460 
470 
480 
490 


KS  = 
KI)  = 
KD  = 
YG  = 
AVOL 
QR  = 
XI  = 
SC  ” 
XC  •» 
REM 
QA  = 
QB  =•■ 
QC  ~ 
QD  = 

525  REM  . 
PI , P4 , P7„ 

526  REM  . 
LENGTH  OF 

527  REM 

540  P2  * 
550 
560 
570 
580 
585 
590 
600 
6 1 0 
620 
630 
640 


60: REM  ...HALF  REACTION  CONSTANT 
.  06  :  REM  .  .  ENDOGENOUS  DECAY  COEFF 
KD  /  24:  REM  . ..(1/HR.) 


(1/DAY) 


.6:  REM  .. 
=  2870000 ! 
.25:  REM  . 
100  :  REM 
50  : REM  . . 
0 :  E  M  .  . 


YEILD  COEFF. 


. SET  RECIRCULATION 
...IN  IT.  COND.  ORG.  CONC. 

I N I T I AL  SUBSTRATE  ( MG / L ) 

I  NIT. COND. 

..FLOWS,  (GPM) 

<Q3  /  1440)  *  1000000! 

(QB  /  1440)  *■  1000000! 

(Q5  /  1440)  *  1000000! 

(QD  /  1440)  *  lOOoOOO! 

.  .P2,,  P5,P8„  PIARE  THE  SPACING  BETWEEN  PULSES. 
PO  ARE  WHERE 

..THE  PULSES  BEGIN.  P3, P6„ P9, AND  PJ  ARE  THE 

THE  PULSE. 

.  ANOXIC  TOWER: 


P3  =  2.4 
P4  =  2! 

P5  »  2 , 5 

REM  TRICKLING  FILTER 


2.  4 


645 


PS 
P9 
F‘0  “•  2 
PI  =  .5 
PJ  ~  2.4 
H  =  .  1 


650  GOSLJB  1000: 
660  GO SUB  2000 : 
665  IF  X$  =  "Y" 

670  GOSLJB  3000: 

671  GOTO  790 
680  GO SUB  4000 
790  LPRIN  r: LPRINT 
810  LPRINT 

820  LPRINT 

821  LPRIN! 

825  LPRINT 

826  LPRINT 
830  I  .PR 


TOWER 

FILTER 


REM  ANOX  1C 

REM  . TRICK. 

THEM  680 

REM  .  .  .  AC  T .  BLlJ.OGi 


L..PR  I  N  T :  L  PR  I  NT 

TAB ( 15) : "TO  CHANGE  hHY  C( JIM. . 


i  AMI'S  OR  INITIAL 


900  END 
1000  REM 


TAB  15)  ;  11  CONDI  'HOPS,,  LI ! 
f  AB  •;  1 5 )  5  "  T  O  CHANGE  Pi  ILSF 
. ANOXIC  TOWER 


at'  ! . I.  NEC  I  0-) 

CIS  !  t.  j.NEi 


G  : 
61 
T  = 

T  A 


1 005 
1 007 
1 0 1 0 
1  o  :i  i 
1020  00 

1022  QO 

1023  00(0) 


QB 

I NT (00 
00 


X  00 


5)  / 10( 


H 


64  0 l: 
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1025  NS  =  SN 
1030  FS  =  SMI 

1032  LPRINT 

1033  LPRINT  TAB (15) 5"  ANOXIC  TOWER" 

1034  LPRINT 

1035  LPRINT  TAB (1 5) ; "TIME  N-INFL.  N-EFFL.  FLOW" 

1037  LPRINT  TAB  (15)  ;  "  (HRS.  )  (mq/L)  (mq/L)  ( GPM ) 

1038  LPRINT  TAB  (15)  5  " - - - - . . . - . - 


1039 

LPRINT 

1 040 

LPRINT  TAB ( 17 

) ; 

T ; TAB ( 28 ) 5  NS $ T AB < 3 / ) 5 FS \ TAB ( 46 ) 5 

1047 

MS 

=  1 

1 050 

FO 

R  I  =  1  TO 

481 

1052 

T 

=  T  +  H 

1053 

tfo 

*  T6  1 

L054 

IF 

T  >  LP  THE 

N 

1450 

1 055 

GO 

SUB  1500 

1060 

K 1 

-  H  *  DM 

1070 

FS 

pq  +  a  g 

•ft' 

K 1 

1080 

GO 

SUB  1 500 

1 090 

K2 

=  hi  *  DM 

1 1 00 

FS 

=  FS  •+•  .  5 

*1 

2  -■  .5  *K1 

1 1 1 0 

GO 

BUB  1500 

1 130 

K3 

=  H*DN 

1 1 40 

FS 

=  FS  +  K3 

... 

5  *  K2 

1150 

GOSUB  1500 

L 1 60 

K4 

~  H  *  DIM 

1165 

FS 

=  FS  -  K3 

i.  170 

FS 

=  F  S  +  ( 1 

/ 

6  )  *-  (  K1  2  *  ( |<2  +  103)  +  K4) 

1. 206 

IF 

MB  =  5  THE 

N 

1 208 

1207  GOTO  1230 

1208  MB  =  0 


1209 

NS  =  I NT (NS  *  100  +  . 

5) 

/ 

1 00 ::  F 

S  -■  I  NT  ( 

FS 

1 00  -i 

7  i  oc: 

) 

1 2 1 0 

00  =  INT(  QO  *  100 

4 

„ 

5)  /  1 

00 

1211 

T  --  I  NT  (  T  *  100  +  . 

K-j 

) 

/  J.  00 

1215 

00(16)  “  QO 

1217 

GOTO  1225 

1 2 1 9 

LPR I NT : LPR I NT I LPR I NT 

1220 

LPRINT  "TIME" ; TAB ( 10) 

n  II 

N- 

I  NFL.  " 

>1  1  hB  (  2U  ) 

5  "M 

... 

EFF'L,. 

" ; TAB (31 ) 5 "FLOW" 

1 2it’'5 

LPRINT  TAB (17) ; T ; TAB ( 

28 

)  p 

NS  5  T AB 

(37)  ;;  FS; 

TAB 

( 46 )  ?  QO 

1 227 

IF  T  >  L.P  THEN  1450 

1 230 

MS  =  MS  +  1 

1 250 

NEXT  I 

1 450 

RETURN 

1500 

REM 

1520 

IF  T  >=  (Pi  +  G  *  P2 

G 

*  F'  3) 

THEN  154 

o 

1530 

GO  1 U  1560 

1 540 

IF  1  M (PI  +  G  *  P2 

+ 

(G 

+  i ) 

*  F'3 )  THEN 

1580 

1560  MB  =  SN 
6 1 5 7 O  GOTO  1 590 
1580  MS  =  SN  +  SI 

1590  IF  T  >=  (P4  +  G1  *  PS  +  G1  *  P6)  THEN  1640 
1630  GOTO  1660 

1640  IF  T  <  =  (P4  +  Gl  *  P5  +  (01  +  1)  *  P6)  THEN  1680 
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1660 

(30 

-  QB 

1 663 

QO (76)  = 

QO 

1 670 

GOTO  1700 

1680 

QO 

=  QA  + 

QB 

1685 

QO ( 76 )  = 

QO 

1700 

IF 

T  (P4  + 

( G 1  +  1) 

*  P6  +  ( G 1  4  1 ) 

1 720 

1710 

GO 

TO  1730 

1720 

61 

=  G 1  + 

1 

1 730 

U 1 

=  (QO 

/  VOL 

)  *  ( NS 

-  FS) 

1731 

U 1 

=  01  * 

8 .  02 

1 :  REM  . 

. . CONVERSION 

1 732 

02 

=  (  - 

ZLAM 

*  (  z 

F()  )  /  (  (QO  /  A) 

1 740 

DM 

~  111  4 

(  02 

*  FS ) 

1 800 

RE 

TURN 

2000 

RE 

T  TRICK 

L.  I  MG 

FILTER 

200 1 

LP 

RINTSLP 

FLINT: 

LPR  I  NT 

2002 

72 

ss<) 

2005 

6 

=  o 

2006 

G 1 

=  0 

20 1 0 

SF 

=  Cl 

20 1 5 

MS 

=  1 

20 1 6 

74 

=~  0 
r  c 

2030 

QT 

-  QD  + 

QO  ( 7 

4 , 

2035 

QT 

(0)  =  Q 

T 

U ) 


2036  Q'l  .  "r(QT  *  jf.iO  f  .5)  /  100 

2037  LP. ,  -IT  TAB  (28)  ;  "TRICKLING  FILTER" 

°038  LPR I NT : LPR I  NT : LPR I NT 

2039  LPR  I  NT  TAB  ( 1 7 )  5  "  T I  ME  C .!.  NFL .  "  s  7  AB  <  35 )  ;  "  C-EF  FL  „ 

FLOW" 

2043  LPR  I  NT  TAB  ( 15)  5  '•  (FIRS.,  )  (it, g/L)  "  TAB  (34)  5  "  (,ng/L  ) 

(GPM) " 

2047  LPR I NT 


2054 

LF1 

•"<  I  NT 

TAB  ( 

17); 12 

; TAB (27) 5  SO; 

2060 

F  01 

9  .J  == 

1  TO  481 

2065 

T2 

=  T2 

4-  H 

2066 

T4 

=  74 

4-  1 

2067 

IF 

72  > 

L.P 

THEN  2 

300 

2070 

GO SUB  2 

500 

2080 

l<5 

=  H 

*  DC 

2 1 00 

SF 

•••-■  SF 

+  . 

K  h 

21 10 

00: 

::>(..)  B  2 

500 

2120 

K6 

=  H 

¥■  DC 

2130 

SF 

=  SF 

•i  „ 

5  ■*  K.6 

—  „  5  V-  KK;« 

2140 

GO! 

SUB  2 

500 

2 1 50 

K7 

-  H 

#  DC 

2 1 60 

SF 

=  SF 

=  K 

7  .  „  5 

•fr  K  6 

2180 

GO! 

SUB  2500 

2 1 90 

I<8 

=  H 

*  DC 

2195 

FS 

F’  3 

—  j-:;; 

7 

2199 

WF7 

/  :r~  X- 

*  ( 

K6  4  |  „. 

7  } 

2200 

SF 

SF 

+ 

1 6666  7 

*  (K.-j  +  WH  i 

2202 

SF 

(14 ) 

SF 

2205  IF  ITS  *  5  THEN  2215 


GOTO 
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22.10 

GO ! 0  2260 

2215 

M3  =  0 

2220 

IF  SF  >  10  THEM 

2235 

2230 

SF  =  10 

2235 

SF  =  I NT (  SF  * 

100  +  .5) 

)  / 1 00 

2  3  / 

SF"  <T4)  ='  SF 

2240 

T2  =  I NT (  T2  * 

1 00  +  . 5 ) 

2250 

L.PR.IMT  TAB  ( 17)  5 

T2  5  TAB ( 27  > 

2260 

M8  =  MB  •+  1 

2270 

NEXT  J 

2300 

RETURN 

2500 

REM 

2520 

IF  7 2  >=  <P7  + 

G  *  F‘8  +  c- 

2530 

GOTO  2560 

2540 

IF  12  <“~  (F'7  4- 

G  *  P8  +  ( 

2560 

SO  =•  CF 

2  5  /  U 

GOTO  2590 

2580 

SO  =  CF  4-  C 1 

2590 

IF  T2  >=  (P7  + 

(  G  +  1  > 

2595 

GOTO  2600 

2599 

G  =  G  +  1 

2600 

IF  T2  >=  (PC)  + 

Gl  *  PI  4- 

2630 

GOTO  2660 

2640 

IF  T2  <*  (PO  4- 

Gl  *  PI  + 

2655 

QT  (T'4)  =  QT  , 

2660 

QT  =  QD  +  Q0(T4) 

2680 

QT  ~  QD  +  QC  + 

QO ( T4 > 

2685 

QT  ( 'T'4 )  as'  QT 

2700 

IF  "l"2  >—  (PO  + 

(  G  1  +  1  )  •;» 

2720 
27 1 0 

GOTO  2730 

2720 

Gl  as  Gl  +  1 

2730 

!J 3  =  (  QT  /  CVOL)  -4  (SO  - 

2740 

U4  =  Cl... AM  * 

(ZT  B)  > 

2750 

DC  as  i.J3  +  (U4  * 

SF) 

2800 

RETURN 

3000 

REM  . . . ACT 1MATE 

1)  SL.U.D(:jfc 

30 1 0 

I  3  —  u 

30 1 5 

MS  as  ;i. 

30 1 6 

QE  Q  T  ( 0 )  4-  <  0 

R  »  QT ( 0 ) ) 

30 1 7 

QE  -  I NT (  QE  * 

1  0  0 ..  5 )  /  1  (')  t 

3020 

S  as  S  C 

302 1 

LPE I  NT 

-.>022 

L PR I NT 

3025 

X  =  xc 

3026 

LPR1NT  TAB < 15) ; 

II 

’■?’  (-)  V '  7 

L.  PRINT 

30.28 

L.. PRINT  FAB  ( 15)  ; 

" T I ME  C 

FLOW 

3029 

LPE 1 NT  TAB (15) 5 

" (HRS. ) 

( OFT! ) 

M 

3030 

LPE I NT 

TAB  ( .1 

c:j  »  it . . . . . . . 

3031 

j. PRINT  TAB  (17)  5 

!  3  |i  1  AB  (  2.6  ) 

/  IQOsQT  =  I  NT  (•  -QT  *  100  4-  .5 

/  J.  00 

)  :!  b'Cl  ■!  T  AB  ( 3 6 )  5  SF-  II  !  AB  ( 4-6 )  |!  Q  I 

3  *  P9)  THEN  2540 

(G  +  '  1)  *  P9>  THEN  2580 

*  P9  +  (  G+l  J  *  P8)  THEN  2599 

o j  *  P J >  THEM  2640 

(Gi  +  1)  *  PJ)  THEN  2680 

PJ  4  (  Gl  -K  j. )  *  PI)  THEN 

SF)  *  8.021 
/  (  (QT  /  CA  >  •  W  > 


AlJ  !  .1.  v'A  !  E.'u  SLUDGE" 
EG .  CONG ERR..  „  CONG  „ 
(mg/L)  >, inq  ■'!....) 

. . .  . . _ . . mtmm . . . . . . . .  i  i 

H  X  n  «  !“J.H  (  O’  /  )  5  "O  |S  I  AEl>  (  4t)  )  'I 
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3050 

FOR  V  =  1  70  481 

3052 

T3  =  T3  +  H 

3053 

T5  =  T  5  •+■  1 

3054 

IF  7  3  >  L...P  THEN  3400 

3055 

S8  =  S 

3060 

X8  =  X 

3065 

I  =  1 

3070 

GOSUB  3500 

3075 

38  =  S  +  S  <  I ) 

3080 

X8  =  X  +  X  (I) 

/  2 

3085 

a  2 

3090 

GOSUB  3500 

3095 

SE!  ~  S  +  S  < I ) 

/  £ 

3 1 00 

X  8  =  X  +  X  ( I ) 

/  ■■  ~t 

/  x'.. 

3 1 1 0 

1  =  3 

3 1 20 

0  0  S  U  B  3 5 0 u 

3 1 30 

SB  S  4-  S  (  1  j 

/  2 

3 1 40 

X8  X  +  X  (I) 

J  A 

/  • 

/ 

3 1 60 

|  .... 

GOSUB  3500 

3170 

8  =:  fe  4  (  1  !  ti ) 

*  1 5  ( I ) 

+  8  (  4  )  4- 

2  *  <  S  ( 3 ) 

3  ( 2 )  )  > 

3180 

X  =••  X  4-  ( l  /  6 ) 

*  <  X  ( 1  > 

*•  X  ( 4- )  4- 

2  *  (  X  ( 2  ? 

4  X  ( 3 )  )  ) 

3 1 90 

IF  M8  -----  5  THEN 

32 .1  0 

3200 

3210 

32 1 5 
3240 

9070  3290 

M8  =  0 

IF  S  >  5  THEN 

■  24'  > 

S  -=•  INI  <8  *  1. 00  4-  .5  )/ 

•|  .  ,.  .■”  x  = 

I N  T  (  X  * 

i.  Ou  +  „  5 ) 

3. 00 

3250 

T3  -  INI'  (T3  * 

1  00  4-  . !; 

)  /  .loo; 

PE  -  INK 

QE'  100 

3270 

LFRIMT  TAB (17) 

•!  1  3  |l  1  f  -'lB  (  2 

6  >  ;;  X ;  T  AB 

(37) ; S; TAE 

(46) ; QE 

3290 

MO  =  M8  4-  1 

33 1 0 

NEXT  V 

3400 

RETURN 

3500 

REM 

3506 

3507 

IF  X  <  10000  THEN  3510 
x  =  loooo 

35 1 0 

QE  =  QT ( T5)  4 

(PR  *  PI ( 

•v;:::-  \  •. 

I  } 

3520 

XO  =•  ((XI  *  QT 

(T5) >  +  ( 

X  *  QR  > 

qt  rr  □) ) ) 

/  QE 

3521 

IF  £7  >  5  THEN 

3530 

ju!..::! 

S  =  5 

3530 

5  <  I )  =  H  *  (  (  ( QE  /  fVv'OI . 

)  ( SF  ( 

i  5)  ~S  > ) 

*•  8.  02 1  -- 

( S  / 

(S  +  KS>  )  *  X) 

i 

3538 

05  -  (QE /  A VOL 

)  *  (XO  -• 

X  > 

4  8 »  02!  1 

3  lC'j  3 

U6  ~  (  ( K  •*  YG> 

*  <  S  /  ( 

b  ■  *  K  S )  > 

*  X)  --  (K 

D  *  X ) 

3540  X  < I )  -  H  *  ( U5  +  06) 

3550  RETURN 

4000  REM  ...ACTIVATED  ULUOGE  FROM  f'RlCKL. IMG  FILTER: 
4100  T4  -  0 
4110  MO  =  1 

4  1 20  QE  -  QT  •  O )  4-  <  PR  07  <  0 )  ) 

41.30  PE  ”  I  NT  ( QE  *  1.00  +  .50/100 
41.50  LPR1NT 
41.60  I  ..PR  I  NT 
4170  X  ~  XC 


* 
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4180  LPRINT  TAB (15);”  ACTIVATED  SLUDGE” 

4190  LPRINT 

4200  LPR I  NT  TAB  (15);  "  T I  ME  ORG  „  CONC EFFL .  CONC . 

FLOW” 

42 1 0  LPR  I  NT  TAB  ( 1 5  >  ;  "  ( HRS )  ( mq  /  L )  (mg  /  L ) 

(GPM)  " 

4220  LPRINT 

4230  LPR I NT  TAB ( 1 7 ) ; T4 ; TAB ( 26 ) ; X ; TAB ( 37 ) ; BF ; TAB ( 46 ) ; QE 

4240  FOR  V  =  1  TO  481 

4250  T4  =  T4  +  H 

4 2 6 0  T  6  =  I  &  +  1 

4270  IF  T4  >  LP  THEN  3400 

4280  S8  “•  SF 

4290  X8  -  X 

4300  I  =  1 

4310  GO SUB  4600 

4320  88  =  SF  +  SF ( I >  /  2 

4330  X8  --  X  +  X  ( I )  /  2 

4340  I  =  2 

4350  GOSUB  4600 

4360  88  =  SF  +  SF ( I )  /  2  - 

4370  X8  =  X  +  X  ( I )  /  2 

4380  I  =  3 

4390  GOSUB  4600 

4400  88  =  SF  +  SF ( I )  /2 

4410  X8  =  X  +  X ( I >  /  2 

4420  I  =  4 

4430  GOSUB  4600 


4440 

SF 

;= 

SF 

4* 

(  .1.  /  A )  ■X'  ( •“4” 

( 1 

1 

•i- 

SF  ( 

4)  + 

2  * 

< 

SF  (3) 

(2)  )  > 

4450 

X  < 

=  X 

+ 

( 1 

/  6 )  *  (  X  (  1  5 

+  X 

( *i- 

) 

+  2 

*  ( 

X  ( 3 

) 

+  X  ( 2 ) 

4460 

IF 

MS 

— 

5 

THEN  4480 

4470 

GO 

TO 

454 

0 

4480 

MS 

as 

0 

4490 

IF 

SF 

5 

THEN  4510 

4500 

SF 

er 

4510 

SF 

I  NT 

f 

SF  *  100  .5 

.) ./  :i. 

00 

1! 

X 

I  NT 

\  X 

* 

1 00  -I- 

4520 

T  4 

=s 

INI 

( T 

4  *  100  +  .5)/ 

1 00 

“ 

QE 

== 

I  NT  ( 

Jjl  !••• 

•I 

l 

00  +• 

4530 

L.F'1 

R I N 

T  T 

AB 

( 1  7 ) ;  !  4 ;  i  AB  ( .26 

t  5  X 

5  T 

AB 

( 37 

)  5  b  F 

5  1  A  B 

( 4 

6)  QE 

4540 

MS 

- 

M8 

4* 

1 

4550 

NE 

XT' 

u 

4560 

RE 

TUR 

N 

4600 

REM 

46 1 0 

IF 

X< 

ioooo  Then  4630 

46-30 

QE 

ss 

Q  !  ( 

TE 

i )  +  ( QR  -ft  Q !  (  i 

5)  ) 

/Q 

E 

4  6  6  b 

XO 

( ( 

X I 

•s-  Q  1  i  5  >  )  +  (  X  * 

Q 

R 

*QT 

•;  T5 ) 

)  )  / 

Q 

E 

4640 

IF 

SF 

5 

THEN  4660 

4650 

SF 

i:4- 

4  A  c;>  (J 

SF 

( i ) 

H 

4  (  'i  QE.  /  A  VOL 

)  * 

( 

SF 

(  T  6 

) 

SF )  > 

* 

8 ..  02 1 

(K* (SF/ 

(S 

F  + 

l< 

:s ) )  *  x>> 

4  c>  /  U 

US 

~ 

( QE 

/ 

A  VOL )  *•  ( XO  ~ 

X) 

* 

8 

»  02 

l 

4680 

LJ6 

= 

(  ( 

K 

*  VG)  * (SE  /  (S 

F  + 

K 

b ) 

)  -ft* 

X) 

--  (K 

D 

X ) 

46  VO 

X  ( 

I  ) 

=  H 

(05  +  U6) 

4700  RETURN 


I  »\J  E>  L.J  5E3  T  Ft  J  f*  L 


W  S  T  EZ  1 R:  ET  A  T  i  -I  E  IM 1  P  8 ^  IM  T 

s . I  o  9 _ E3  T  O  IM  *=%  FT  F:> 


HOW  MANY  HOURS  WOULD  YOU  LIKE  TO  SEE 
8 


FLOW  FROM  AREA  B  (MOD)  IS  4.5 


FLOW  FROM  AREA  A  (MOD)  IS 


ENTER  VALUE  FOR  COD  (rnq/L.)  IN  TRICKLING  FILTER  INFLUENT 
300 

DO  YOU  WANT  THE  IR1CKL I  NO  FILTER  EFFLUENT'  TO  RE  THE 
ACTIVATED  SLUDGE  INFLUENT  ?  ( Y  or  N> 

Y 


ANOXIC  TOWER 


TIME  N-INFL. 
'•HRS,,  )  (mq/t.. ) 


N EFFL.  FLOW 

(mq/L)  (OEM) 


0 

•••  ;.i  r::; 

3 1 25 

i;-. 

••■"i  tz: 

.c*.  •_.! 

1 .,  38 

3125 

1 

2! 

J.  „  07 

•v  ^  •;>  l-j 

1 ,.  5 

~r 

1 . 07 
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TRICKLING  FILTER 
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APPENDIX  D 

NPDES  PERMIT,  1980,  SUMMARY 
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